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We construct the path integral for one-dimensional non-linear sigma models, starting 
from a given Hamiltonian operator and states in a Hilbert space. By explicit evaluation 
of the discretized propagators and vertices we find the correct Feynman rules which differ 
from those often assumed. These rules, which we previously derived in bosonic systems [lj, 
are now extended to fermionic systems. We then generalize the work of Alvarez-Gaume and 

- i— i i 

Witten Q by developing a framework to compute anomalies of an n-dimensional quantum 
field theory by evaluating perturbatively a corresponding quantum mechanical path integral. 
Finally, we apply this formalism to various chiral and trace anomalies, and solve a series 
of technical problems: (i) the correct treatment of Majorana fermions in path integrals 
with coherent states (the methods of fermion doubling and fermion halving yield equivalent 
results when used in applications to anomalies), (ii) a complete path integral treatment of 
the ghost sector of chiral Yang-Mills anomalies, (Hi) a complete path integral treatment of 
trace anomalies, (iv) the supersymmetric extension of the Van Vleck determinant, and (v) 
a derivation of the spin-| Jacobian of Alvarez-Gaume and Witten for Lorentz anomalies. 
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1 Introduction 



In this article we discuss Euclidean path integrals for one- dimensional systems with 
a Hamiltonian which is more general than H(p, x) = T(p) + and their applica- 

tions to anomalies. Namely we shall consider non-linear sigma models whose classical 
action is of the formf] J « ^9ij^f^-dt, plus fermionic extensions, which may be, but 
need not be, supersymmetric. This is the area of quantum mechanical path integrals 
in curved space, a difficult and controversial subject [§, |J. 

Quantum mechanical path integrals are of importance as toy models for path 
integrals for field theories. Because the former are finite, no renormalization is neces- 
sary and subtle issues can better be studied. In addition, quantum mechanical path 
integrals are useful because some quantities of n-dimensional quantum field theories 
can be calculated in a much simpler way by using the corresponding one-dimensional 
path integrals. The prime example are anomalies, which we shall discuss in detail 
in the second part of this article. In the first part we give a careful derivation of 
quantum mechanical path integrals and the Feynman rules to which they give rise. 

This article is an extension of a previous article Q on bosonic quantum mechanical 
path integrals to the fermionic case and to applications to anomalies. In section 2.1 
we briefly review the bosonic path integrals, but most emphasis in section 2 and 
beyond is on the fermionic case. At the end of this introduction we shall state which 
results are new, but we shall start with a general introduction in which we discuss 
bosonic and fermionic systems on equal footing. 

The basic problem we solve is the following. Given a Hamiltonian H(p, x, ■0) 
with arbitrary but a priori fixed ordering of the bosonic operators pi, x l (i = 1 . . . n) 
and fermionic operators^ ip a , ip\ (a = 1 . . . n), find a path integral representation for 
the transition element (z, f}\ exp (^—j-H^J \y, x) where \y, x) is an eigenket of x % and ip a 
with eigenvalues y % and x a , while (z, fj\ is an eigenbra of x l and with eigenvalues 
z % and 77 a . Following Dirac || and Feynman M we shall insert N — 1 complete sets of 
x-eigenfunctions, N complete sets of p eigenstates and N complete sets of coherent 
states 

J \x)y/rtx)(x\<Px = J \p)(p\d n p = J \r,)e-™ a (f)\ f[(dfj a dv a ) = 1 (1) 

a=l 

6 Hats denote operators, but where no confusion arises we will omit them. 

7 Curved indices in space-time will be denoted by /i, v, . . .. In the non- linear sigma model the 
corresponding indices will be denoted by i,j, . . . 

8 For Majorana fermions, we shall formulate one approach where a — 1, . . .n and another where 
a = 1, . . . , n/2. We shall find it convenient to use fermionic operators satisfying anti-commutation 
relations without h, {ip a , ip},} = $b ■> w bilc = ihSj as usual. 
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to obtain a phase space path integral in the limit iV — > oo of the form 

(z,f]\exp (^-J-H^J \v>X) ~ J dx i dp i dfj a dr ] a e~^-e Ldt (2) 
where L = —ipji^ + hfj a f] a + x, fj, rj). However several questions arise: 

(i) which is the relation between the operators H(p, x, iji) and the functions 
H(p, x, fj, 77)? Different operator orderings of H must lead to different functions 
H. Are there particular orderings of H for which H is equal to the naive classi- 
cal Hamiltonian? After integrating out the momenta, is the action in the path 
integral invariant under the usual Einstein (general co-ordinate) transforma- 
tions or supersymmetry transformations if the corresponding Hamiltonian H 
commutes with the generators of these symmetries? (The answer is no). 

(ii) What are the Feynman rules needed to evaluate the path integral perturba- 
tively? 

(iii) What is the precise meaning of the measure dx l dpidfj a dri a ! Is there a normal- 
ization constant in front of the path integral, or even factors g(z) a g(yY where 
g = det(g i j) f ! (The answer is yes) 

(iv) By adding couplings to external sources, one obtains propagators. These propa- 
gators are not the usual translationally invariant propagators because they must 
satisfy boundary conditions at t — —f3 and t — 0. What boundary conditions 
for Pi(t) and the fermions must be imposed? 

(v) When one is dealing with Majorana fermions, how should one define coherent 
states, and what is the Hilbert space in which H is supposed to act. Should one 
impose boundary conditions both at t — —f3 and t = for Majorana fermions, 
and if so, how is this compatible with the fact that the equations of motion are 
linear in time derivatives? 

These are some preliminary questions. Not all of them are new, but we shall 
automatically get answers by following our derivation of path integrals. These an- 
swers will be summarized in the conclusions. The most important question we solve 
has to do with Feynman rules. The propagators (x l (a)x j (r)) and (ip a {<t)^1(t)) for 
configuration space path integrals (and in addition {x l (a)pj(r)) and {pi(&)pj(j)) for 
phase space path integrals) with — 1 < er, r < contain factors 9 (a — r), where a and 
t are the time divided by (5. In configuration space, there are contractions of x l (a) 
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with ar 3 '(r), which contain a factor S(a — t). When one computes Feynman graphs, 
one has to evaluate integrals over products of these distributions, for example 

/ = J° J 8{a - r)6{a - t)6{t - a)dadr (3) 

In addition, there are in general equal-time contractions with factors 6(0) and 0(0). 
How should one evaluate such integrals? The function 8(<j — t)6(t — <j) is nonvanishing 
only at a = r, a set of zero measure, so any smoothing of 8(o~ — r) while keeping the 
product of the two #'s intact leads directly to zero. A more natural way would seem 
to be to expand 5(a — r) and 8(a — r) into an infinite series, and to cut off these 
series at some large N ('mode regularization'). (The series expansion of 8(a — r) 
could for example be defined so that d a 6(a — r) = 5(a — r)). Performing the whole 
calculation for finite N, one would expect to obtain the correct answer by taking the 
limit N — > oo at the end of the calculation. This is incorrect as we shall see but 
first we must answer a fundamental question: what does the expression 'the correct 
answer' mean? 

Several authors have tried to give meaning to continuum path integrals in curved 
space, in particular configuration space integrals of the form /[cfe 1 ] exp(— hS), by 
freely inventing definitions which maintain Einstein invariance at intermediate steps. 
The transition element (z,fj\exp (—jrH^j \y,x) being in principle known from either 
heat kernel methods [0, |^| or direct operator approaches^ the loop calculation based 
on these path integrals should in the end reproduce the results for the transition 
element. Here problems arise: no prescription for L is known which keeps covariance 
at all stages and which gives the correct result at two- or higher loops. 

Our point of view is the following: we define the continuum limit path integral 
as the limit of the discretized path integral obtained from (|]). Hence, for us 'the 
correct answer' means: the answer which reproduces the results of the Hamiltonian 
approach. Integrals over products of discretized propagators, vertices and equal-time 
contractions are well defined and finite and by taking the continuum limit, the correct 
Feynman rules automatically emerge. The results are that 5 (a — r) should still be 
considered as a Kronecker delta function, even in the continuum case. So, for example, 
the correct value of I in (H) is 1/4. To bring out the surprising consequences of these 
new Feynman rules, consider another integral 

I = J° J° 5(a - r)9(a - r)9{a - r)dadr (4) 

9 Writing (z\ exp (—^H^J \y) as J (z\ exp(—j-H\p)(p\y)d n p, it is clear that by expanding the expo- 
nent and moving all pi to the right and all x % to the left, keeping track of commutators, coefficients 
of a given power of (3 are finite and unambiguous. Similarly when fermions are present. 
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If one were to require that 5 (a — r) = d a 6(a — r) even at the regularized level, one 
would find 1/3 for this integral, whereas the correct answer (and the answer obtained 
from treating 5(cx — r) as a Kronecker delta) is 1/4. This demonstrates that mode 
regularization is incorrect. 

No expressions with higher powers of 5 (a — r) arise if one introduces in config- 
uration space path integrals the extra ghosts of ||. These ghosts are necessary to 
exponentiate the factors g 1 ^ 2 which arise if one integrates out the momenta. These 
factors g 1 ! 2 were first found by Lee and Yang in a study of non-linear deformations 
of the harmonic oscillator |T0j and were written by them as extra terms in the action 



of the form g 1 / 2 5(0). We find it much more convenient for higher loop calculations 
to replace them by local terms with ghosts, similarly to the familiar Faddeev-Popov 
ghosts in gauge theories. The phase space path integrals are always finite since the 
propagators do not have any 5 (a — r) singularities, and with these ghosts also the 
configuration space path integrals are finite^]. 

In the first part of this article we shall give a careful derivation of the new Feynman 
rules. We begin by rewriting the Hamiltonian H in Weyl-ordered form (H)w, and 
then use that matrix elements of exp(— jr(H) w ) with e = (3/N can be immediately 
evaluated to order e by using the 'midpoint rule'. (Weyl ordering is discussed in 
section 2.3.) Namely one may replace exp(— |if) by exp(— f (H)w) in the kernels of 
the path integral 

J (a; fc ,7/ fc |exp(-^if)|p fe ,£ fc _i)e - ^ (5) 

and then one may replace in {H) w the operator pi by (pk)i, x l by \{x\ +x l k _ 1 ), ipl by 
&-i,a, and r by + Vi-i)- (Or $ by i(jj fci0 + £fc-i,o) and ^ by The net 

effect is that one can extract the function H{p kl ~(xk+Xk-i)i £fc-i> 5 l+^fe— i))- F° r 
linear sigma models (with H = T + V) rigorous proofs based on Banach spaces and 
the 'Trotter formula' |TT]] exist These do not apply to non-linear sigma models, 
but we have found a simple derivation of (§) which is precise enough for our taste. 
Note that no matter whether H is gauge invariant or not, this midpoint rule holds; 
it is a purely algebraic result. 

By using the background field formulation and coupling quantum deviations to 
external sources and decomposing the action S in a suitable free part and an 
interaction part, we find discretized propagators and vertices in closed form. The 



10 Power counting would seem to indicate that there are linear divergences due to the double 
derivative interactions. This would seem to contradict the theorem that quantum mechanics is 
finite. The ghosts save the theorem. 



bosonic discretized propagators were already found in 0, while fermionic propagators 
for Dirac fermions were already found in W2 . 



We shall first consider Dirac fermions, but then an even number of Majorana 
(real) fermions ip a (t) and operators i[j a satisfying the Dirac brackets {tp a , tp b } = S ab . 
To define coherent states we need creation and annihilation operators, and these 
we shall construct in two different ways: by doubling the number of fermions by 
adding a second set ip% of free fermions, or by 'halving the number of fermions' and 
constructing ip A and ipA from pairs il) 2a ~ l and ip 2a as (ip 2 ^ 1 ± iip 2a ) / \^2. In either 
case we Weyl order, use the fermionic midpoint rule, and find propagators. The 
propagators (ip a (a)ip k '(r)) are different depending on whether one doubles or halves 
the Majorana fermions, and also the action in the path integral and the transition 
elements differ. In applications to anomalies, however, these differences disappear. 
In the conclusions we explain the reason for these results. 

In all cases (bosonic systems, Dirac fermions, Majorana fermions either with 
doubling or halving), the propagators can be factored as follows 

P 



[z, rj\ exp 




~ 1/2 

where Sdass is the classical action, D s contains the Van Vleck determinant and takes 
care of the one-loop contributions, while to order j3 A n contains the trace anomaly 
for n = 2 dimensions (it is proportional to the scalar curvature R). This structure of 
the propagator was proven for general bosonic systems in 0. For fermionic systems 
Ds is actually the superdeterminant of — -^rSd^s-^j where $ J = {z l ,fj a } and $j = 
{y\Xa}- The fact that the one-loop contributions should be equal to the Van Vleck 
superdeterminant is a check on the correctness of our Feynman rules. 

In general, the action in the path integral contains extra terms of order h and h 2 . 
These extra terms are due to rewriting the Hamiltonian in Weyl ordered form. In 
our case, we shall only encounter ft 2 terms. Hints that such terms might be necessary 
in the path integral were first found by DeWitt J7|. Schwinger, who studied the 
Poincare operator algebra for Yang-Mills theory in the Coulomb gauge (which is a 
non-linear sigma model in 4 dimensions), found that one had to add extra non-naive 
terms of higher order in H to the Hamiltonian in order that the algebra closes ||13||. 



Subsequently many others have studied these extra terms JT4[ [TIJ. To check our new 
Feynman rules and also check that no further modifications at order h 3 are present, 
we perform a 3-loop calculation in appendix A.l. 

Our methods also apply to systems with more than two momenta (higher deriva- 
tive theories). We consider in appendix A. 2 such a system, and check that the results 
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of the phase space approach agree with those of the configuration space approach. 
This is Matthews' theorem It provides another test on the correctness of our 
results. 



The second part of this article contains applications of these quantum mechanical 
path integrals to anomalies. Anomalies of an n-dimensional field theory can according 



to Fujikawa be written as [|17 



An = Tr(je~P R ) (7) 

where J is the Jacobian for a symmetry transformation of the path integral, and R is 
a regulator. A general method to construct consistent regulators R which maintain a 
given set of symmetries is given in |18[]. It was first proposed by Alvarez-Gaume and 



Witten to consider a corresponding linear or non-linear sigma model in one dimen- 
sion, for which R becomes the Hamiltonian H. The basic idea is that the Fujikawa 
trace can be viewed as the trace over the Hilbert space of a quantum mechanical sys- 
tem with a finite number of operators x^jd/dx^, Dirac matrices 7° and, if present, 
internal symmetry generators T a . After representing the Dirac matrices and internal 
symmetry generators (by means of Majorana fermions and an auxiliary ghost system) 
as operators in the same Hilbert space, the Fujikawa trace can be rewritten in terms 
of a suitable quantum mechanical path integral. The path integral can be used to 
compute matrix elements, and in terms of those the anomaly becomes 

A n = J dx i ^Jg(x)dr] a df] a e flaVa {x,fj\Je~K H \x,r}) (8) 

and by inserting a complete set of states between J and exp(—/3H/h) one obtains a 
product of the transition element and the matrix element of J. The anomaly is the /3- 
independent term. Depending on the anomaly, i.e., depending on the matrix element 
of J, there are different factors of (3 in front of this expression, and since (3 counts the 
number of loops, different anomalies require a different number of world-line loops to 
be evaluated. 

The simplest anomalies are the chiral anomalies. For these (|8]) is actually (3 
independent, due to the topological nature of chiral anomalies, and, as we shall show, 
as a result one only needs to evaluate one- loop or tree graphs. In fact, Alvarez-Gaume 
and Witten wrote the whole expression in (|8|) again as a path integral of the same kind 
as we consider for the transition element, but now with periodic boundary conditions 
both for the bosons and for the fermions. The one-loop contribution for such path 
integrals can then easily be written as the determinant of the kinetic operator of 
deviations about classical solutions. We shall obtain, of course, the same results but 
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will start from the transition element we obtained previously, and then simply do the 
rest of the integrals in (|S|). In our approach the ghosts for internal symmetries are 
part of the complete path integral, and are not treated by operator methods and by 
projecting on one-particle states as in 0. The results of Alvarez-Gaume and Witten 
were extended to trace anomalies in || . In this case the product of the Jacobian and 
the transition element could again be written as a path integral with now anti-periodic 
boundary conditions for the fermions, but since the background fermions are constant 
if they satisfy the equations of motion, the authors of || had problems in finding the 
correct boundary conditions at t = —(3 and t = for Majorana fermions. Rather, 
they used an operator formalism for the fermions, and operator-valued actions. We 
shall present a complete path integral formulation. Since the trace anomaly receives 
contributions from higher loop graphs, the details of the path integral do matter very 
much. For example, forgetting the extra h 2 terms in the action or the measure factor, 
one obtains incorrect results. We shall also give a derivation of trace anomalies. 
Although we will not do so here, one can use the same framework to derive the 
gravitational anomalies due to spin- 1/2 and spin-3/2 fields. The only subtlety is 
the question precisely which Jacobian one should use in ([?]) for spin-3/2, and this is 
discussed in detail in appendix A. 3. 

We conclude this introduction by stating which of our results are new. The 
Feynman rules for products of distributions^ are new; for the bosonic case they 
were obtained in |]| and for the fermionic case here. Our treatment of coherent 
states for fermions (first considered in [^C]) follows [21] and is somewhat simpler (we 
believe) than those treatments in the literature which use four kinds of coherent states 
(namely bras and kets which are eigenstates of $ or tp^). For a good discussion of the 
latter see p2| . The careful treatment of Majorana fermions (doubling and halving), in 
particular the fact that the transition elements are different, is new, as is the proof that 
in all cases the one- loop contributions sum up to a superdeterminant. New in section 
three is the complete path integral treatment of chiral and trace anomalies, with 
no need to introduce matrix-valued Hamiltonians or to perform certain projections 
on the ghost states by hand. Also new is the complete diagrammatic evaluation of 
these anomalies in section 3 and appendix A. 4. The final new result is the correct 
incorporation of 'Lee- Yang' ghosts for higher derivative theories in appendix A. 2. 

Besides all these new results, we have spent a great deal of time to make the whole 



11 Mathematically, it is possible to multiply objects called generalized functions that contain the 
set of distributions. However, to determine which generalized function corresponds to a given dis- 
tribution still requires an 'underlying physical principle' (see e.g. [jl9| and references therein), and 
our rules can be seen as an example of such a principle. 



subject of quantum mechanical path integrals for non-linear sigma models accessible 
to a large audience, and we hope that this article will also be a useful review of known 
results. 

2 Path integrals for finite time 

In this section we establish the framework we need in order to perform calculations 
using one-dimensional path integrals for finite time. Given a quantum mechanical 
system defined by a Hamiltonian with a certain ordering prescription, we derive the 
corresponding path integral formulation. This includes both the action to be used 
in the path integral, as well as the Feynman rules, the latter being not only a set of 
expressions for the propagators and vertices, but also the correct prescription how to 
evaluate integrals over products of these, as they occur in actual loop calculations. 
We will first discuss the bosonic non-linear sigma model, and derive the correct rules 
for configuration space path integrals. For phase space path integrals for bosonic 
systems, see [0. Then we will derive similar results for the extension to complex 
(Dirac) fermions, for which we will need to introduce coherent states. Finally we 
shall discuss the modifications which one must make for Majorana fermions. Along 
the way, we will present a variety of checks that our Feynman rules are the correct 
ones, by evaluating various transition elements and comparing these with the results 
obtained from operator methods. Further convincing evidence is provided in appendix 
A.l and A.2. 

2.1 Bosonic non-linear sigma model 

We will start by computing the transition element for the following quantum mechan- 
ical Hamiltonian 

H = \g- 1/A PiV9g l3 Pig- 1/A (9) 

This operator in Einstein invariant if Pi is hermitian. (The p t transform under Einstein 
(=general co-ordinate) transformations as p\ = f {f^;Pj}> if the inner product is 
defined by (|l|) , from which the Einstein invariance of H follows ) . The more general 
case, where also a scalar and vector potential are present, can easily be found at 
each stage in the computation by covariantizing the expressions, and including the 
scalar potential in the interactions. The cases of other Hamiltonians, for example 
Hamiltonians whose operator ordering is different from (^), will also be clear. 
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One can evaluate the corresponding transition element, 

T(z,y;/3) = (z\ew\-jrHj \y), (10) 

for finite (Euclidean) time f3 through a direct, but rather tedious, computation in the 
operator formalism p3| , p], |24]. Namely, by writing T as / d n p(z\ exp(— $H)\p)(p\y) 



and expanding the exponent, moving all pi to the right and all x % to the left, keeping 
track of all terms with up to two commutators, we find the following result (correct 
through order /5, counting z — y as being of order /3 1//2 ) 

T = (2vr^)-"/ 2 exp {-\s^P]) D 1 ' 2 exp ("^) (H) 

where (2whP)~ n / 2 is the usual Feynman factor, 5^ [2, y; /5] is the classical action for a 
geodesic with x{—0) = y, x(0) = z, — j^R yields the trace anomaly in n = 2 dimen- 
sions^ and D 1 / 2 is proportional to the square root of the Van Vleck determinant |p5| 
which gives the one-loop corrections to the transition element 

D 1 ' 2 = P^g-^z) det ^^S4z, y] P]^ 2 g-^(y) 

= l-^(z)(z-yY(z-yy + 0(p 3 / 2 ) (12) 



The results in (|TT|) and ( |l2l) agree with DeWitt's classic paper in which he uses 
heat kernel methods. Note that if one views T for (3 — > as the kernel of a continuum 
path integral action, then D 1 ! 2 corresponds to a non-local term. Each of the factors 
in ( |TID is a general co-ordinate bi-scalar (a scalar both in y and in z) . We stress that 
the answer for T in flTID is finite and unambiguous. 

We will now construct the path integral whose loop expansion reproduces ([TT|). 
First we rewrite the Hamiltonian in Weyl-ordered form ||26|| , which for any monomial 
in p and x is defined by [n + m)\(p n x m )w = dad™(ap + bx) m+n [^7|, ^8], ^9f. For the 
Hamiltonian in (0) we find the well known result ]27L |3U| 

# = (\g*PiPj) w + y (^ r K* + R) (is) 



Then we use the correspondence between Weyl-ordering and the midpoint rule [31. 27 



d n p(x k+1 \G w \p)(p\x k ) = d n pG(p,x k+1/2 )(x k+1 \p)(p\x k ) (14) 



12 Our convention for the curvatures are R{T) pcrp v = d p Y +T pT v T ap T — (p <-> cr) = R{lu) p<J a b e^e-a 
with i? (9cra 6 = dpUJaab + ^ pa c ^crcb - (p <-> c) and = i?(r) A1CT „ CT , so = R(uj) ppa be ap , while 

i? = . 9 ^i?^ = (<^v CT - gwsrfadvgpa + .... 
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where we denned Xk+1/2 = fak + ^fc+i)/2, and G is any function of p, x. Clearly, 
exp (-f# w ) = (exp (-f#)V + ^(e 2 )^ so usin S (0) we find 

y cf> exp f-^fl^ |p){p|»*) 

= J d n p exp (-|fl"(p, ^+1/2)) b) (pkfc) + 0(e 2 ) (15) 



Note that i?(p, Xk+1/2) hi (jig) contains the terms in (p~3| ) of order ft . We argued in 
U that the terms denoted by C(e 2 ) will not contribute to the path integral, and can 
therefore be neglected. This is the only point in our derivation of the path integral 
that is not mathematically completely rigorous. What is subtle is the meaning of 
0(e 2 ). For example, one can view p either as being of order 1 or of order e -1 / 2 
(since there are Gaussian integrals with exp(— p 2 /2e)). In the latter case the terms 
denoted by 0(e 2 ) are actually of order e. Each case leads to a different kernel. We 



argued in |I] , using the effective potential trick of fllj] and ]32| , that both kernels are 
equivalent under the path integral. However, this is not a rigorous argument. Ideally 
one should keep all terms on the right hand side of fll5]) that might contribute in the 
limit iV — > 00, evaluate the path integral at the discretized level and then prove that 
in the limit iV — > 00 all extra terms do indeed drop out. 

We now insert N — 1 sets of x-eigenstates and N sets of p-eigenstates into 
(z\ exp — £-H\y), and we arrive at the discretized phase space path integral using 



Jd n x J g(x)\x)(x\ = 1 = Jd n p\p)(p\ and (z\p) = (27rfr)~ n / 2 (exp jrp ■ z)g~ 1 / 4 (z). In- 
tegrating out the N momenta we find the discretized configuration space path in- 
tegral, with iV factors g 1 ' 2 (xk+i/2) i n the measure from the p integrals, iV products 
g~ 1 l i (xk+i)g~ l l 4: {xk) from the inner products (x\p) and N — 1 factors g 1 ^ 2 (xk) from 
the completeness relation in x-space. The action is given by 



N-l 



fc=0 



1 h 2 e 



—9ij(xk+i/2)(xk+i - x k ) l (x k+1 - x k ) J + — (TV + R){x k+ i /2 ) 



(16) 



where we define x^ = z and x = y, and e = (3/N. We decompose x 3 k into a 
background and a quantum part, and S into a free and interacting part 



bs.fc 



+ ql S = S (0) + S (int) ; k = l,...,N-l (17) 



where 5 (0) = J2 k =o j- e gij(z)(q k+ i - q k ) l (q k+ i - q k y ■ We take the metric in 5 (0) at 
z in order to facilitate comparison with (0), although any other choice should give 
the same result. (Of course, propagators and vertices will be different if we make 
a different decomposition into a kinetic and interaction part, and also the measure 



factor (see (22)) will be different, but the final results should not change). Since 
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we take Xb g to be a solution of the classical equations of motion of S^, in general 
g(mt) con t a i ns j n addition to the true interactions also a pure background piece and 
terms linear in q 3 k . However, counting z — y as being of order /3 1//2 , we need only 
a finite number of tree graphs and tadpoles at a given order in {3. The N factors 
<? (2^+1/2) are exponentiated following || (for an alternative approach see [33]) by 
using anti-commuting ghosts b and c and a commuting ghost a 



exp 



" 2f3 2 h aij ( Xk+1 / 2 ^ (^+1/2*4+1/2 + a k+l/2 a k+l/2) 



Since the constant K will cancel, we do not determine it, and the reason for the 
particular normalization of the ghost action will become clear later. Introducing 
modes for the quantum fluctuations q by the orthonormal transformation 

^E^^in^); k,m=l,...,N-l (19) 

we may change dx k — *> dq k — > dd? m . Obviously, the Jacobian for this transformation 
is 1. The quantum part of the action becomes equal to 

1 N ~ 1 1 Tfl 7T 

- J-S (0) (q) = - E - cos—). (20) 

Next, we couple to external sources 

^source) = _ £ £ Qk+l ~ Qk + Q^^^ + SQUrces for ghosts \ (2 1) 

k=0 ^ € ' 

so that we can extract the exact discretized propagators in the usual way. Completing 
squares and performing the final integration over d 3 m , b° k+1 , 2 , <4 +1 / 2 , a J k+1 / 2 leads to 
N — 1 factors g~ l l 2 (z) and N factors g x l 2 {z) as well as an overall factor (27rfr/?)~ n / 2 . 
(The factor (2ixht)- Nn ' 2 which comes from the p integrations and the normalization 
of the plane waves combines with the factor Ilm=i( 7re ^) n//2 (1 — cos-^) - ™/ 2 from the 
Gaussian integrations over d^ m to yield {2 r Kh(3)~ n ^ 2 since Ilm=i2(l — cos 2 ^") = N). 
Hence 

T=(^^j (2vr^)- n/2 exp(-^ (int) )exp(-^ (source) ). (22) 

The measure factor (g(z)/g{y)) 1 ' i is due to the split S = S<® + S (int) where in S<® 
the metric is taken at z, gij(z) (cf. question (iii) in the introduction). Different 
splits clearly lead to different measures, but we continue with (p2"D. Such factors 
are often ignored but are crucial to obtain the correct transition element. Equation 
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( ^2]) is to be read as usual in path integral formulations, namely S^- mt ' contains only 
derivatives with respect to the sources and background fields, while 1 5'( source ) is the 
function bilinear in sources that appears after doing all the integrals, and in the final 
result we are supposed to put all the sources equal to zero. Thus, S^ nt ^ contains 
the discretized vertices while 5 ,< - source ) yields the discretized propagators. Defining 
ik+1/2 = (xk+i — Xk)/e and omitting superscripts and a factor of hg^(z) for the time 
being, the propagators come out as follows 



< <Zfc+l/2ft+l/2 > 

< <?fc+l/2<?Z+l/2 > 

< <?fc+l/2<?«+l/2 > 

< 6fc+l/2Q+l/2 > 

< a fc+l/2<^+l/2 > 



-^(2* + 1)(2J + 1) + |(2nun(fc, I) + 1 - 



+ 9, 



k + 1/2 
JV 

"JVC + ~e 5h > 1 

Ok,l 



k,l 



(23) 



where O^ i is a discretization of the 9 function: 9k,i — if k < I, 9kj = 1/2 if k — I and 
9 kji = 1 if k > I. 

As an example, let us give a more detailed derivation of the < ^+1/2^+1/2 > 
propagator. The partition function (we set the sources F equal to zero, and suppress 
internal indices) reads 



N-i e 

z i g ] = / n ^ ex p~i" 
i=i n 



N-l 



fc=0 



1 ^g fc+ i - g fc V 



G 



fe+1/2 



/ gfc+i + ?fc 
V 2 



(24) 



Now make a change of variables as defined in fllTf ) and (20), complete the squares 
and do the Gaussian integrals over the modes d m . Up to an overall numerical factor 
which was already taken care of in (|22|), the result equals 

c3 N-l 1 



Z[G] 



exp 



— Y - 

2Nh f-J 1 



3=1 



— cos(jn/N) 



/N-l 



2 -, 



X; G fc+1/2 sin((A; + l/2)j7r/JV) cos(j7r/(2JV) 



,fc=0 



(25) 



Differentiating with respect to the sources G gives rise to the following expression for 
the propagator 



< <?fc+i/2<?m/2 > 



2iV £ |_ 



COS" 



; (j7r/(2JV) 



sin((fc + 1/2)j'tt/JV) 



sui(j7r/(2JV)) 
sin((Z + l/2)j7r/^] 



sin 



(W(2JV) 



(26) 
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Each term can be written as a sum of powers of e ln ^ 2N \ and performing the sum over 
j yields the propagator as given in (|2lf). The remaining discretized propagators can 
be found in a similar fashion. We require that x 3 hgk satisfies the boundary conditions 
and the equation of motion of S^ '. In the continuum limit this becomes x J hg (t) = 
zi + [z — y)H/(3, while cp(t) vanishes at the endpoints. In this limit the two-point 
functions become (reinstating the superscripts, factors of g^(z), and defining t = (3t) 

<q i (a)qi(r)> = -{3hg^ (z)A(a, r) 

<b\a)c i (T)> = -2(3hg ij (z)d?A{a,T) 

<a\a)a>{T)> = phg* (z)<%A(*,t) 

A(a,r) = <t(t + 1)6{<t-t)+t{(j + 1)6(t-(j). (27) 

Note that A(er, r) = A F (a - r) + or + \(o + r), where A F {o -r) = \{o - t)9(ct - 
t) + |(r — <j)9(t — a) is the Feynman propagator, and formally d%A(a, r) = 5 (a — r) 
while A (a, r) = at the boundaries. However, the 5 (a — r) is a Kronecker delta and 
moreover the equal-time contractions^ are unambiguously defined. Kronecker delta 
here means that J dx8(x)f(x) = /(0), even when / contains a product of 9 functions. 
From (|23|) we further find in the continuum limit 

< q\^W{r) > = -f3hg lj {z){a + 0(r - a)) 

<?W(t)> = -(3hg^(z)(l-5(a-r)). (28) 

All propagators are now proportional to f3h (this motivated the normalization of the 
ghost action in (fL8|)), and the interactions are given by 



}_ c(int) _ _{_ f° 

h pnJ-i 



^gij(x hg + q) {(x bg + q) l {xhg + 4Y + 6V + aV} 

+ph [ u hrr + R)dr - ls(°\ 

J-i 8 n 



dr 



dr (29) 



(3h J-i 

Clearly, the interactions only depend on the combination f3h~. 

To compute the configuration space path integral, we note that we must ex- 
pand the measure factor g 1 ' i {z)/g 1 ' 4 (y) in (|22| ) and evaluate all vacuum graphs with 
external Xb g , using the propagators in (0), (|28"D and the vertices in (p9|). The q- 
independent part of S( mt > does not yield the full of (|TT| ) since Xb g is only a so- 
lution of the equation of motion; rather, tree graphs with two vertices from 

13 Equal-time contractions in quantum field theory can in general only be fixed by imposing a 
symmetry principle |34[ |. In our case they are fixed by our requirement that the path integral 
reproduces the Hamiltonian results. 
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g(int) con tribute to order (3 terms of the form ^(dg) 2 (z — y) 4 , see fl66l). In the one- 
loop graphs with one vertex S^ 111 *) one finds equal-time contractions proportional 
to (z — y) k dkgij times < q l qi + b % c? + a l a^ > in which the 6(0) cancel, yielding 
J ^^ J° 1 a(d a d T A + d 2 A)dadr = — ~, which cancels a similar contribution from the 
non-trivial measure factor. There are many other one-loop and two-loop graphs, and 
the contribution of each corresponds to a particular term in ([11]). In particular, the 
two-loop graph with one qq, one qq and one qq propagator agrees with (|TTD only if one 
uses /_j J° 1 6 (a — t)6(o — t)8(t — a) = 4, in agreement with the discretized expressions 
for the propagators in fl23|). Adding all contributions we have found complete agree- 
ment. The non-covariant vertices ^(TT + R) conspire with the non-covariant vertices 
found by expanding gij(x) and yield the Einstein invariant expression ([TT|) . The Feyn- 
man rules one has to use in this calculation follow from (^), and they amount to 
the following. First, one writes down expressions for all Feynman diagrams using the 
propagators given by (^) and (|28|) . Adding everything, all divergences coming from 
products of delta functions will cancel (the ghosts of || are crucial for this). The 
resulting integrals should be worked out using the rules that delta functions should 
really be seen as Kronecker deltas and that 6(0) = 1/2. If there are explicit delta 
functions in the integrals, one should be careful with partial integrations and identi- 
ties like /* /' = f(b) — f(a), since these are not always compatible with our Kronecker 
delta prescription PJ. Luckily, in practice we never need to partially integrate. 



2.2 The fermionic case 

We now repeat the analysis of the last section for the fermionic case. We will work in 
a basis of coherent states, and the derivation of the path integral is analogous to the 
one in phase space for the case of the bosonic non-linear sigma model. We consider 
operators i[) a ,ipl, a = 1 . . . n, satisfying the anticommutation relations {ij) a ,$1} = 5%. 
These operators ip are obtained from the canonical variables by rescaling with a factor 
offt~ 1/2 . As a consequence, terms of the form f^Rip 4 or hcuip 2 are terms of the classical 
action, not higher-loop terms. For cases such as the N=l supersymmetric non-linear 
sigma model, where only Majorana fermions are present, we will need to replace these 
Majorana fermions by Dirac fermions. This will be discussed later, here we will derive 
the general expression for the path integral with Dirac fermions. 

Coherent states are defined by 

\v) = e#"\0) ; (i?| = <0|e** (30) 

satisfying ^77) = r]\rj) and (77 = (fj\fj. We could in addition also define coherent 
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states build around a Dirac vacuum (completely-filled Fermi sea) (rj\ = 
(— l) n (0|^ n • • • ip 1 (e^^j , satisfying {r]^ 0, = (^|^ a , and similarly \ff) = 
(e^ fe )^i...^t|0)(-) n . These states are often used in the literature for the con- 
struction of fermionic path integrals in a way which closely mimicks the x ■ p ap- 
proach, see e.g. |22fl . In contrast, in our approach we only use the coherent states 
(|30l); although both approaches are completely equivalent, we believe ours is more 
economical. The inner product and decomposition of unity read formally the same 
as for bosonic coherent states [|2l] 



<tj|0=e* ; 1 = dfjd£\$e-*(r}\ (31) 



but note the ordering of the anticommuting variables. Our convention is that dfj = 
dfj n . . . dfji, while d£ = d£i...d£ n , or equivalently, df/dC, = Y\k=i(dVkd£,k)- Hence 
J <i££„ . . . £1 = 1, and J d£ II5b=i £fc = (— l) n//2 for even n. With these conventions, the 
trace of an operator over the fermionic Fock space is given by 

trace (A) = J dtdfje* (fj\A\£) (32) 

Again, we define Weyl-ordering by (n + m)\(ip n ip^ m )w = d^d^(f]ip + r]^p J< ) m+n , with 
the fermionic derivatives acting from the left. For an arbitrary Weyl-ordered operator 
G we can now derive the midpoint identity for coherent states 

(v\G\v) = JdxdxG(x,^Y^-)e~ xx (v\x)(x\v) 

dxdx G{^—^-, X)e~ xx (v\x) (x\v) 

d X d X G(XiX + (1 - \i)fj, A2X + (1 - h)v)e-™(v\x) (x\v), (33) 

with A1A2 = 1/2. This formula can be proven in the following way (cf. [0). It is easy 
to check that it is valid for an operator G of the type (ip J< ) k , which is automatically 
Weyl-ordered. Then one uses that, for a Weyl-ordered operator A, the operator 
B = (tpA ± Aip)/ 2 (the sign depending on whether A is bosonic or fermionic) is 
also Weyl-ordered, and hence that any Weyl-ordered operator can be obtained by 
repeatedly applying this identity to an operator of the type (ip^) k . We can now 
proceed by inserting unity in (|33|) for G = B, using its validity for A as induction 
hypothesis. We find 

(fj\B\ V ) = Jdxd X e-^^(fi\x}(mv) 

d x d X d^e-^^(f}\ X )A^,^)e-^(mm 

d^B(l^L)e-tt(v\0(£\v) (34) 
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where we have used that J d\dx e~* x J '(x) = /(O)- This completes the proof of the 
first two lines of (|3"3|). The third line of ( |3"3"D can be demonstrated in a similar fashion. 
We can now apply this identity to G = exp (— \Hw) > an d use that G is Weyl-ordered 
to order 0(e 2 ). We find, neglecting these higher order terms, after inserting unity 
N - 1 times, with f3 = Ne, 

(fj\exp(~H w j\x) = 



N-l 

n 

fc=i 



d£kdi k e exp ( -t^w ) 16) 



ft 



(6| exp |x) (35) 



where we defined £ 



(£ fe+ i | exp ( -j&w 



rj. Now use the midpoint rule for each matrix element 

16) = 



di/j k dipk e-^^+5 fc+ i^+^a. exp I ~-H{ij k 



(36) 



We can integrate out the £,£ i n ( HI ) (first the £, then the £) to obtain 

P 



exp 



ft' 



-if 



N-l 

Y[ dip k dt/} k exp 

k=0 



Tjlp 



N-l 



N-l 



fc=o V 



T lpk-lpk-l , Iq-zt ^h + ^k-1- 
e ft V 2 



(37) 



where we defined = x- The first term in the integrand is the usual boundary term 
one obtains in path integrals for coherent states (cf. |21], 0). From this result we 
conclude that the action in the continuum path integral is J®p(h$ip+H)dt—hijj(0)ip(0) 
with the boundary conditions V>(0) = fj and 0) = X- Notice that the boundary 
term is essential to produce the correct equations of motion x/j — ip — 0. After 
decomposing ip k and i(j k into a background piece and a quantum piece, we couple the 
latter to external sources. Putting all of H into H mt , we obtain the discretized v/ji/j 
propagator by inverting the kinetic term matrix Aj )k = 5j yk — 5j,k+i- The result reads 



-1 k < I 
k > I 



(38) 



If we now define i> k -i/2 = (ipk ~ V'fe-i)/ 6 an d ^k-i/2 = (V'fc + ^fe-i)/2, we obtain 

< i>ki*i-i/2 > = 



< i> k ipi-i/ 2 > 



-Oi,k 

i 



(39) 



where we recall that 9 k ,i = if k < I, 9 k ,i — 1/2 if k — I and 9 k< i — 1 if k > I. We 
can now, as in the bosonic case, write down the corresponding continuum expressions 
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which we will use in actual computations. However, it is important to realize that in 
diagrams in which products of S and 9 functions arise, we are now able to resolve the 
resulting ambiguities by returning to the discretized expressions (|39|). In particular, 
we find again the 'rules' that 9(a,a) = ~, and Jdadr9(a — t)9(t — a)5(cr — r) = |. 
The continuum propagators read 

< $ a {a)ip b {r) > = -5 ab 9(r - a) 

< $ a {a)ip b {r) > = -5 ab 5(a - r) (40) 

There is no factor of h because we chose to work with operators ijj a , ipb satisfying 
{ip a , ipb} = 5%. Clearly, in this derivation we could instead equally well have started 
from the second line in (|3~3"D and have introduced ipk+1/2 = (rfk+i ~ V ; fc)/e, 4>k+i/2 = 
+ V'fc)/2, leading again to (^+1/2^) = — @i,k f° r the discretized propagators, 
and therefore also to identical Feynman rules in the continuum limit. 



2.3 Weyl-ordering of N=2 and N=l Hamiltonians 

We will now derive the Weyl-ordered Hamiltonians corresponding to the supersym- 
metric N=2 and N=l Hamiltonians. These are the most interesting fermionic sys- 
tems, and play the same privileged role as (||) in the bosonic case. As it turns out, 
different expressions result when we take the independent fermionic fields to have flat 
or curved indices, and Weyl-order with respect to these independent fields. Let us 
first consider the N=2 quantum Hamiltonian 

n = Pi-^iabVa^i (41) 

where is the spin connection and a = 1,2. We can now define Dirac spinors in 
the following way 

r = ^w+^2) ; r = ^«-# 2 a ) (42) 

satisfying the anticommutation relations {ip a , ip b } = S ab . In terms of those, the N=2 
quantum Hamiltonian corresponds to the following classical action 

C = \g l3 (x)x l x 3 + hr(r + iVVfc) - hi 2 Rabc^W^ h ^ d - (43) 

Its field equations read 

-j^-C = -9ijjf j + %R(u) ijab x j ^ b -h 2 D l R(oo) abcd riJ b iJ c ^ d 



Sip a dt 

() C = ^ a -fiR(uj) dbc a ^ d ^ c , (44) 



5ip a dt 



where 



and 



—jj = x j + T j kl x k x l (45) 

(JjL 

Ey* = ip" + xWb^, = ^ + i'ui^t. (46) 

The invariance of the action under the two rigid supersymmetries follows by contract- 
ing the field equations with the variations 

6x* = eeif + fele 
SiJ; a = x i e'*e-5x i u i ab i; b 

Stfj a = -eefx* - 5x i u i a %. (47) 

All terms then cancel (using the cyclic and the Bianchi identities for DiR abc d). Since 
this action is in Euclidean time, it is not hermitian, nor is (Sifj a Y equal to 5ip a , but 
it can be obtained from a hermitian action in Minkowski space by the Wick rota- 
tion tu = —He- The classical Noether charge for supersymmetry reads eia(x)ifj^,x l 
with a = 1,2, while the quantum charge is given by Q a = e^x^^^nig^ 1 ^ = 
g~ 1 l i Tt i g 1 l i e ia (x)il)a- ft is hermitian and Einstein invariant, and {Q a ,Qp} = 25 a pH. 
This shows that the Hamiltonian is supersymmetric and that supersymmetry is pre- 
served at the quantum level. Note that the variation 5ip a + 5x l Ui ip b is covariant; 
the 'pull-back' term 8x l Ui ip b is due to the presence of fermionic equation of motion 
terms in the bosonic equation of motion. 

Weyl-ordering of the bosonic part of the Hamiltonian gives rise to the usual 
contribution fllBD 

^ 2 (r + ^"rjrjj (48) 

For the Weyl-ordering with respect to the fermions, we first choose the fermions with 
flat indices as our independent variables. Clearly the terms quadratic in the fermions 
yield no contribution, because their anticommutator is proportional to S ab , which 
gives zero upon contraction with Ui ab . For the terms quartic in the fermions, we use 
the identity (see e.g. |j36| ) 



i { [r, , [r, } = (^v W) w + \^ hc (49) 
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Since the four fermion terms in the Hamiltonian have exactly the same symmetry 
as the operator on the left hand side of the above identity, we can easily deduce the 
fermionic contribution to the Weyl-ordered Hamiltonian to be 

-ltf(R + gij U J> U ./) (50) 

Adding this term to the bosonic contribution (^8[) we find, for the N=2 case, 

\h 2 <ft (r£r} fc - u ia b u jb a ) (si) 

Alternatively, one can take the fermions with curved indices, namely ip % = ^(V'i+m/'I) 
and ipi = gij-j^(ip{ — itpi), as independent variables, and Weyl-order with respect to 
these. One finds then that now the bosonic and fermionic contributions exactly cancel, 
or, in other words, the N=2 supersymmetric Hamiltonian expressed in these variables 



is already Weyl-ordered |36|| . 



The N=2 Hamiltonian cannot be interpreted as the regulator of a corresponding 
quantum field theory, because each of the 2n ip* (a — 1, 2) would have to correspond 
to a Dirac matrix, whereas there are only n Dirac matrices in an n-dimensional 
quantum field theory (with x % with i — 1 . . . n). However, it plays a role in the path 
integral evaluation of the index of the <9-operator of the Dolbeault complex || . 

We now consider the N=l supersymmetric case, where only one species of Ma- 
jorana fermions is present, which makes the generalization of the previous result 
non-straightforward. The Hamiltonian in this case is equal to 



H N=1 = y^g^jg-^-'tfR 



= Pi- —uj iab ip a ip b (52) 

where a — 1 . . . n, and the fermions satisfy the usual relation {ip a , ip b } = 5 ab . 

This Hamiltonian cannot be obtained by truncation of the N=2 Hamiltonian. 



For example, putting ipi — ip 2 = requires e\ + 62 = 0, see (f47|), but the resulting 
Hamiltonian has —j^h 2 R instead of —^h 2 R, and is no longer supersymmetric. The 
reason is that the truncation ipi — 1^2 = is no longer consistent at the quantum level 
since {ipi — tp 2 , ipi — ^2} is nonzero. The easiest way to obtain (|52|), is to start from 
the N—l action with £ = ^giji; l x^ + \^ a %^ a -, to construct the Noether quantum 
supersymmetry charge Q = g 1 ^e l a ^p a ii i g~ 1 ^ 4: , with Hi = Pi — Y UJ iab' l P' :i "4 )b ', and then to 
evaluate H = \{Q, Q}. The algebra is the same as used to evaluate {Ip , Ij) } in ( |1 1 1| ) , 
and leads to (^2|) . 
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In order to construct coherent states, we cannot work with Dirac brackets or 
Majorana spinors, but we need creation and annihilation operators. There are two 
ways to achieve this: either by combining interacting Majorana spinors or by adding 
free Majorana spinors. 

We will first combine these Majorana fermions into complex spinors x> X i n the 
following way 

y A = ±W"- x + ; x A = -j=W> 2A - 1 -i1> 2A ) (53) 

where A — 1 . . . n/2 and {x A , X B } = 8 AB - The inverse relations are given by 

i) a = ^(x (a+1)/2 + X (a+1)/2 ) if a odd 
v2 

r = -±=( x ^-x a/2 ) if a even (54) 

We substitute ( |54"D into ( |52"D and define Weyl-ordering with respect to the operators 
X, X i n the usual way. It is easy to prove, considering separately the cases a, b odd or 
even, that the following equality holds 

r^ b = (r^) w +\s ab (55) 

We can now derive the Weyl-ordered expression corresponding to the N=l Hamil- 
tonian. The bosonic part yields the same contribution as before, see (f48|), and, 
because of the above identity, the part quadratic in the fermions is again already 
Weyl-ordered. It remains to consider the term quartic in the fermions ff(q uartlc ) = 
— ^T? g % i uJiabU) j C< i r i\) a '^) h, ^) c '^) d . If only one pair of fermions gives rise to a non-trivial anti- 
commutator, we obtain a contribution proportional to g l iLdi a c LVj C bi/j a ij) b , which vanishes 
identically for symmetry reasons. We therefore need also the second pair of fermions 
to yield a non-vanishing anticommutator in order to find a non-zero contribution. 
Now there are two possibilities: the two sets of fermions are in different sectors (have 
different Dirac index A), in which case we only need to Weyl-order both pairs sepa- 
rately using (p3), leading to 



+ -g ij u ia b u jb a (56) 

where we should remember that, if we define A = [(a + l)/2] etc., not all four indices 
A, B, C, D are identical. The second possibility corresponds to the case A, B, C, D 
all identical. In that case a, b, c, d are all equal to 2k + 1 or 2k + 2, and we only need 
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to consider the case that two of them are equal to 2k + 1 and the other two equal to 
2k + 2. In this case the Weyl ordered expression vanishes^. Hence this leads to the 
same result, so that in fact ([56]) is valid for all a, b, c, d. For instance ip 1 ^ 2 "^ 2 ^ 1 = | 
and (%p 1 %p 2 "ip 2 %p 1 )w = 0; we end up with a coefficient \ in (|56|) because we obtain the 
same contribution from a term ip 1 ^ 2 ^ 1 ^ 2 — — j. Adding the bosonic and fermionic 
parts, we find the total contribution from Weyl-ordering to the scalar potential 

l -h 2 (r + ^rirj fc ) - ^ftVW V" (57) 

As one might have anticipated, the Weyl-ordering of the Majorana fermions yields 
half the result for Dirac fermions, see (|51f) . 

The other way to construct Dirac fermions from Majorana fermions is to add a 
second set of free Majorana fermions. Denoting the original fermions ip a by ip®, and 
the new ones by ip%, we again construct Dirac fermions x an d X ( as i n ®))> but then 
we use the N=2 formulation given before. The four fermion term in the Hamiltonian 
now reads 

h 2 

~ -^t^UabUJjcd^ (58) 
o 

and we should in this case Weyl order it with respect to x a an d X a - Again, (|55|) holds 
for ip%, even though the definition of x a an d X a is now different. Using the operator 
identity 

ip^iplipf = ^(6 ab 5 cd - 5 ac 5 bd + 5 ad 5 bc ) 

+ (^ ab (i'ii , i)w + five more terms) 

+ (<VtoiV (59) 

we find that the two-fermion terms vanish due to anti-symmetry while the double 



contractions yield the same answer as in (|57|) . 



Finally we can also Weyl-order the N=l Hamiltonian with respect to the fermionic 
variables x an d X i n ©)' but now with curved rather than flat indices. In that case 
one expects to be left with a remainder 



: rr(i?+^T*r;. fe ). (en) 



1 To see this, take k = 0. Any real linear combination of ip 1 an d can be written as aip + aijj 
for some complex a. The Weyl ordered expression of an arbitrary combination of ip 1 and ip 2 is 
proportional to the sum over all graded permutations of a, j3, 7, S of the operator (aip + aip)(/3ip + 
Pfy^ip + "f^)(8ip + Sip). This equals a/3jSipi/j + afijS'ipip, and the sum over graded permutations of 
the ordinary constants a, /3, 7, S and a, /3, 7, S clearly vanishes. 
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2.4 Evaluation of supersymmetric transition elements 



2.4.1 The N=2 case 



We will first compute the transition element for the N=2 case using the path integral 
formulation. We shall obtain all terms through order (3. We rescale r = t/(3 to 
make the /3-dependence more explicit and facilitate keeping track of the order in the 
expansion in (3. First note that we can write 



S 



q I okin I nint 
^bos ~T ^->f er T O fcr 



(61) 



Here S^ os contains all terms with only bosonic or ghost fields. Note however that it is 
not identical to the action we wrote for the bosonic path integral, as the Weyl-ordering 
of the fermionic terms also gives a contribution ~ (R + u 2 ). The contributions to the 
path integral involving only this part of the action can trivially be found from the 
bosonic case, because the latter terms can, through order /3, simply be replaced by 
their classical expectation values. The other two terms are given by 

o 

VcWA> fe - W o (0)/(0) (62) 



nkin 
°fer 



and, after integrating out the momenta pi, 

o 

fdr 

-1 



Qtnt 
°fer 



(63) 



In the background field approach we decompose i/)(t) = iphgir) + Vv( r )' and ^(r) = 
"0b g ('7") + V'qu('7")- We choose , 0bg('r) and T/»bg(r) to be the solutions to the equations of 
motion of the kinetic part of the action, (|52]) , which satisfy the boundary conditions, 
i.e. we take ^§ g (T) = ff and V , b g ( r ) = X a - Since we require ip a (0) = ff and tp a (— 1) = 
X a , this implies that the quantum fields need to satisfy the boundary conditions 
$J U (0) = and ^ u (-l) = 0. Then simplifies to j\ % u % u dr - f] a X a , i.e., there 



are no terms linear in fermionic quantum fields in (p2[) . Of course, (p3|) does contain 
such terms. 

We can now compute the transition element by expanding exp (—^S^j and 
contracting the quantum fields, using 



h 



oint 



bos 



(29) + hp j\-R + gVuMVj 1 *) 



(64) 



and the propagators for the fermions given in (^0|). When we expand exp (— 

we will for the first term in this expansion only need the contraction at equal time 
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< ip a {T)%l) h {T) >= —7;h5 ab . Hence, the first term in (|63| ) will yield no contribution, 
whereas the second term contributes \(3hR ab fj a x h — \(3TiR. The term — ^[3hR cancels 
the R term in (pi]). Next consider the term | (hSj^j , where only terms from the 
square of the xuipip term contribute to this order. When we contract one pair of 
fermions, we find the uu term in fl67l). When we contract only the four fermionic 
fields, or two x fields and two fermionic fields, one finds zero. Finally, we can contract 
four fermionic and two bosonic fields. This yields a contribution — ^(3hg^uj ia b ujj b a , 
which cancels the ujuj term in (|64[) . Contractions involving the other bosonic fields 
are again of higher order and need not be considered. 

Taking all contributions into account, we find for the amplitude 



\Z,rj\ exp 



\y, X ) = (2vr/%)- n/2 exp(- 



■—Sb — ttSf 
n n 



1 - —(3hR{z) - —R l3 {z){y - z)\y - z) j + -(3hR ab (z)fj a X l 



12' 



12 



(65) 



where 
Sb 



12(3 



—gij{z){y - z)\y - z) j + —d kgij (z)(y - z)\y - z) j (y - z) k 
d k d Wij (z) - lg mn (^(z)T n kl (z))(y - z)\y - z)\y - z)\y - z) 



(66) 



is the expansion through order f3 of the length of the geodesic joining z and y (cf. 
(O)), and 



S f 



-U ah rfx h ~ h(y - zyu iab (z)v a X b 

—zHv - z T(y - zy(diuj jab (z) + u ta c (z)uj jcb (z)^fj a x 



^(3K'R abcd {z)n a X b ri c X d 



(67) 



The uj and duj terms are obtained by expanding the first term in flB3|). In the contin- 
uum limit, Sf becomes the fermionic action, including the correct boundary term, as 
derived above 



h ldr[ <WAV + x l uo iab r^ b - i-(3hR abc ^ b ^ d ) - M a6 rf<(0). (68) 



We can easily check that the expansion through order (3 of ( |68| ) indeed equals the 
expression in (67) when the equations of motion are imposed. The latter follow from 
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([HP (of course we also need the bosonic part of the action to find the full equations 
of motion) 

■£ + r jk x j x k - ///»",„,>' r'V' 1 + hi 2 g ij [d 3 R abcd ] ^>V C ^ = 
r + x l Ui a ^ b - 7iR\ cd f^ d = 

^ a + x i oj i a b $ b - hR a bcd ^ c '4) d = (69) 

where a dot denotes differentiation with respect to t. We can now expand the La- 
grangian in a Taylor series around its value at t = 0, and then do the trivial time 
integrations. This yields 

S = PL(0) - \(3 2 m + ... (70) 

We thus expand all fields in the Lagrangian around their values at t = 0, making use 
of the equations of motion fl69|). The expansions up to the order we need are given 
by 

^ a (0) = ff 

W (U) - [diujj b - Ui c ujj b ) x 



(71) 



Inserting these expansions into Sf in ( p5| ) and ([T(]) yields the expression (|5T[) . 



We will now show that the final result for the transition amplitude can again be 
written as the product of three factors: a term containing only the scalar curvature 
which is related to the trace anomaly, the exponent of the classical action, and the 
square root of, in this case, the supersymmetric generalization of the Van Vleck 
determinant. The latter we define by 

D s = sdetD AB ; D AB ^ (S B + S F ) ^ (72) 

where $ A = (z l ,fj a ) and <3> B = {y\x b )i an d for S B and Sf we substitute the expres- 
sions (|66|) and (|67D. To evaluate D s write 

Dab = | ^ ** ) (73) 
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We find, expanding in normal co-ordinates around z to simplify the expressions, 



B 



ib 



-gij(z) + hdiU jab (z)r] a x b 

~{diU jah {z) + djU iab (z) + Lu ia c (z)u jcb (z) + uj ja c (z)uJi cb (z)^fj a x b 

-hu iab (z)fj a 



C aj = hu jab (z)x b 



D 



ab 



it 

fi5ab + Kv - z) l ujia b (z) + -(y - z)\y - z) 3 (diUJjab(z) + uJ ia c (z)uJjcb(z)) 



+h(3(R abcd (z) - R, 



adcb[Z))V C X d 



(74) 



We do not need terms of order j3 in B and C, since D$ = det A det _1 (.D — CA~ 1 B) 
and A' 1 is already of order (3. Writing A^ = ^gik(S k j + Pha k j) and D ab = h(5 ab + d ab ), 
we can write the expansion of the super Van Vleck determinant as 



1 + \(3htia + -flhtrCB - hid + ^(trrf) 2 + ]tr(d 2 ) 
2 2 2 8 4 



L Z, Z, 

Multiplying by g~ 1 / 4 (z)g~ 1 / 4 (y) to transform D l J 2 into a bi-scalar, we obtain 



(75) 



Df = m^g-^(z) D y 2 g-^(y) 
1 



1 - -R^(z)(y - z)\y - z) 3 + ^hR abV - a X b 



12 

So indeed we can write 
(z,fj\exp (-j-h) \V,X) = (27rnp)~ n / 2 D 1 s /2 exp 



■{Sb + Sf 



12 



(76) 



(77) 



All terms involving the Ricci curvature in (|65| ) are thus completely accounted for by 
the super Van Vleck determinant, which clearly would not have been the case if we 
had used the ordinary determinant. Finally we note that if we would have rescaled all 
fermions by a factor of (fifi)* 1 ' 2 , then all classical terms are proportional to l/(/3h), 
while all one-loop terms in (|76| ) are f3h independent and the two-loop term in ( [771 ) 
remains proportional to f3H |2j 



2.4.2 The N=l case 



Next we consider the N=l case. We will first evaluate the transition element when we 
double the number of Majorana fermions, and afterwards consider the case that we 
combine the Majorana fermions into half as many Dirac fermions. In the first case we 
add n free fermions ip2 an d combining l/v / 2(V ; i +^2) = an< ^ Vv^V'i — ^2) = 
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we construct the path integral with the corresponding coherent states. The kinetic 
part of the fermionic action is again 

1 °r 

-Sg n = jdr - « a (0)/(0) (78) 

-l 

and yields the free field equations ip a = ip a = 0, but the interaction part containing 
the fermions is now equal to 



i M _ 9 

-S^ = Jdr 

-l 



fj X (79) 



Note that the terms linear in < qu (O) in ([75]) cancel. 

We have introduced an extra set of fermions that do not couple to any of the 
other fields, this way making certain that we do not alter the dynamics. In order to 
preserve local Lorentz invariance we require that the fermions ip% are inert under local 
Lorentz transformations. The extra contribution to the scalar potential from Weyl- 
ordering the N=l Hamiltonian is of the form ^-(IT + R) — jqUJlu to which one should 
add the term —^-R from ([52]) , and its integral can be replaced to order (3 by (3 times 
its classical value. The purely bosonic sector of the path integral can be evaluated 
exactly as before, so we only need to consider the sector involving fermions. We use a 
background field expansion as in the N=2 case, again with constant background fields 
which satisfy the boundary conditions, and substitute tpf = ^(fj a + X a ) + ^l.qu 
([79]) . Since the interactions depend only on ipi, it will be easier to use the propagator 

for Vl,qu 

< ^, q >)< qu (r) >= U a \e{a - r) - 9{r - a)) (80) 

which can trivially be found from the ipip propagator in ([40]) and the propagators 

{ipip) = {ipip) = 0. We now evaluate (exp — hS^). We first consider terms from 

contractions in —-S^*. The equal time contraction of %p1%p\ in ( p0|) vanishes, so this 

term will yield no contribution. Also the contribution from the contraction in x l LJi a b 

o 

in ( |79|) vanishes, as it is proportional to / dr < x J (r)x J (r) > which is zero. Next we 

-l 

consider the term ~ {^S\^ . The contraction of x l (a) with x^(r) is of order j3, and 
would contribute, but its integral over a and r vanishes. One x l (a) contracted with 
a^'(r) is also of order (3, but the other i J would leave a factor x l hg which is of order 
/3 1 / 2 , so this term is of higher order. The contraction of only one pair of fermions 
vanishes for symmetry reasons, so there is now no ujuj in (|8"2|). The contraction of all 
four fermions is nonvanishing, and produces the term jg{z l — y l ){z^ — y : ')uji a b ujji ) a in 



Bl|). Finally, we can contract four fermionic fields and two bosonic fields x % . This 
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yields — j^/3hg^uJi a b u!j b a , and this term exactly cancels a similar noncovariant term 
in the scalar potential due to Weyl-ordering. Adding all contributions we find the 
transition element for the N=l case with fermion doubling 



(z,f}\ex V (-^Hj\y, X ) = {2n^)- n ' 2 exp (~(S B + S F ) 

1 + ^PhR(z) - ^R rj (z)(y - z)\y - 
+ ^(y - z)\y - zyuj ia \z)u ]b a {z) 



where 



— Sp 
n 



s ab v a x b - - A {y - zyu iab (z)(f) a + x a W + xl 

-\{y - z)\y - Z yd i u jab (z)(f] a + X a )(fj b + X h ) 



52) 



and Sb is the bosonic part of the classical action. It is the same result as obtained 
directly from operator methods 



The terms in Sp are obtained by expanding the following classical continuum 
action around z 



Is = -r/> a (0) + f° dt 
n J -i 



1 



J3) 



The equations of motion read 






x l 



1 + r%x j ± k - fJh-x3R^ ab (ip + + 



i 



i) a + ~x i uji ab (ip + ipy 



where Rij a b = diUj ab + u iac ujj cb — {i <-> j). From them one derives further 



x l (0) 



1™ 



(3h 



z- y y- -rUz)(z - y y(z - y f + ?-{z - y yR jab (v + xTiv + xf + 



^(0) = x «-^ z -y) k u k a \r 1 + x) h + \{z-y)Xz-yyd l u j a \r 1 + x) b + ... (85) 

Substituting these results into S = —f] a il) a (0) +L(0) — | jj:L(0) + . . ., the contribution 
from — Ijjr^'i/O cancels the terms in —fj a ijj a (0) + tp a (0)ip a (0) , and one indeed arrives 
at (H). 
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We will now show that the expression for the propagator can again be written as 
the product of the super Van Vleck determinant, the exponent of the classical action, 
and a term involving the scalar curvature which, as shown in section 3.3, determines 
the trace anomaly of a spin-| field. Defining the super Van Vleck determinant as in 
the N=2 case, we find from (1821) 



A 



1.1 



C a j 
D ab 



±g tJ (z) + ~(^u; ia6 (z) - d 3 u iab {z)) (rf + X a ){ff + X b ) 
\u jab {z){f? + X b ) 

h h 

ftSab + 77(2/ - z) l uj iab (z) + -(y - z) l {y - z) 3 diUJ jab (z). 



Using again fl75|) , one finds 



D 



1/2 



1 - ii^O)(y - zY(y - *y + -^(y - - zyu ia \z)u jb a {z) 



where only the last term in (|75|) did contribute and yields the last term in 
we can indeed write 



(86) 

(87) 
So 



;> 7j|exp (-jrHj \y, X ) = {2^y n/2 D 1 J 2 exp {~{S B + ft 



1 + —phR 

24 H 



similarly to the bosonic and N=2 supersymmetric case. 

We will now repeat the analysis for the N=l case when we do not introduce an 
extra set of fermions, but instead combine the n Majorana fermions ip a into n/2 Dirac 
fermions ^ A as ^ A = ^(^ 2A_1 + # 2A ), V A = ^ (ip 2A ^ - i^ 2A ). The kinetic 
part of the fermionic action is now equal to 

1 °f 

-Sg n = J dr6 AB ^ B - 6 AB * A (0)* B (0) (89) 



The interaction part containing the fermions is still equal to 



dr 



^ Qint 
^fer 



'-x i u ia bi/j a i/j b 



(90) 



but now the ip a should be expressed in terms of ty A and ^/ A . We again make a 
background field decomposition as ^/ A = x A + ^q U ; = f] A + ^? A U - Again it will be 
convenient to rewrite the propagators for the Dirac fermions in terms of the Majorana 
fermions. We now find 

1 



< i) a (a)^\r) >= -5 ab (6(a - r) - 6{r - a)) + K 



ab 



(91) 
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where K ab = jd a+1 ' b (l — (— l) a ) — (a <-» 6) (in other words, X is the matrix — |r 2 in 
the 2x2 subspaces). Furthermore, we define the 'background' Majorana fermions 

r = ^( x ( a + 1 )/ 2 + r^+ 1 )/ 2 ), a odd 
v2 

4> a = — ^(x a/2 -V a/2 ), a even (92) 

v 2 

We start again by considering contractions in —hS}^. The equal time contraction of 
ip a ip b is now equal to K ab . These terms therefore contribute 



l -( y - z y Uiab K ab + \ 

2 

qmt < 

2 \h 

we obtain the tree graph 



(y - zfuj iah K ab + -{y - z)\y - zyd t uj, ab K ab (93) 
Next we consider contractions in | \ \S % ^\ . When we contract one pair of fermions, 



hy - z)\y - zYu tab u Jcd (-AK ac ^ d ) . (94) 



8 

When we contract two pairs of fermions, we find the one-loop graph 

- z)\y - z) j u; iab u jcd (^5 bc 5 ad - 2K ac K bd + K ab K cd ^j (95) 

The last term is a product of two one-loop graphs, but should of course not be counted 
as a new two-loop graph. When we contract the two x fields, we only find a contri- 
bution if in addition we contract all fermionic fields. This yields — j^phg^Ui a b u!j b a , 
which cancels again a similar contribution from Weyl-ordering (see fl57D). Adding all 
terms, we find for the transition element for the N=l case without fermion doubling 

(z, fj\ exp (~j:HJ \y, X ) = (2tt^)-"/ 2 E exp {-\(Sb + S F )) [l + ±/3hR] (96) 

where Sp is equal to the background part of the fermionic action (together with the 
boundary term) plus the tree graph of (|94D 

±S F = -6 AB fj\ B -±(y-zyu iab (z)i; a iP b 

-~(y-z)%y-zyd i u; jab (z)ri> b 

+\(y - z)\y - zyoo iab uJ lcd K ac ^ d (97) 

The one but last term is due to expanding u(z + (z — y)r) around z. As the notation 
indicates, to order (3 S F is equal to the classical action (|59"D and (|90| ) with the equations 
of motion satisfied. Their solution is, however, quite a bit more complicated than it 
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was in the previous cases, due to the different boundary conditions we have to impose 
here. All one-loop contributions are contained in E, 

1 



E 



1 



iiRij ( z )(y- z y(y- z y 



+ 



-(y- z) l uj iab + ^(y~ z) l (y - zydiUjab ) K 



ab 



+ -(y- z)\y - zyuj iab uj jcd { -5 bc 5 



bczad 



+ K ab K cd - 2K ac K 



(98) 



Comparing with fl82|) , fl37|) and (^8|) which yield the transition element with fermion 
doubling we note three differences: (i) the background value of the ip in the interac- 
tions are defined differently, namely x a + V a = V^ipf in (0) and ip a in (|2|), (iz) the 
boundary term fjx contains half as many terms in (|97|) as in (0), and (Hi) there are 
extra terms in (|97|) and (p8|) proportional to if and i^i^. Yet, as we shall see, these 
different transition elements yield the same anomalies. 

Motivated by our earlier results, we will now compare the expression E to the 
square root of the super Van Vleck determinant. Defining the super Van Vleck de- 
terminant as before, we find using (|9~7D 

-gij(z) + -(diU jab (z) - d j u iab (z)^ijj a ijj b 

+ ]ji^iab^jcd + UJ jab UJi cd )K aC tjj b ^ d 



BiB 

C A j 
Dab 



B 



1 d 



Sab + 



Ay- z) l u iab (z) + -(y - z) l (y - z) 3 diuj jab (z) 



--uj iac uj jM (z - y)\z - y)°K G 



drj 



I* m - 



dx 



B 



(99) 



The expansion of the super Van Vleck determinant through order (3 can again be 
written as 

,1/2 



9 



l '\z)g^\y) 



l + ~/3tra + -/3trC£ 
i tr( i + i(trd) 2 + itr(d 2 ) 



Now using identities such as 



zAB 







drj 



1 1 '■ 



a i) b 



d 



B 



-2K 



a I) 



(100) 



(101) 
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and 

-— — r f^> 6 ) — — ^7 tt^tt -7— r = 2(R ad R bc - R ac R bd ) + - (5 ad 5 hc - <J ac <? M ) (102) 
ar\ A v ' ax <9?7 v 1 o\ 2 

we find that the super Van Vleck determinant indeed equals E. Hence the same 
factorization which we found in the case of fermion doubling in (|88D also holds in the 

~ 1/2 

case of fermion halving, even though the separate factors D s : and Sp are different. 



3 Anomaly calculations 

In this section we compute the chiral anomalies due to spin-1/2 fermions coupled to 
external gravity and external Yang-Mills fields. Then we consider the trace anomaly. 
In appendix A. 3 we discuss gravitational anomalies for spin-1/2 and spin-3/2 fields. 

In a quantum field theory, the anomalies can be written as the trace over the 
product of a Jacobian and a regulator. Both the Jacobian and the regulator depend 
on which fields one considers as the independent fields, for example, for a scalar field 
4> possible choices are itself and = (fyg 1 ^. Of course, the anomaly itself should not 
depend on this choice of basis. It has become common practice to take <fr for scalars, 
and x a = X ^ 1 ^ 4 f° r spinors as basic variables, because then (with the corresponding 
regulator) the absence of Einstein anomalies becomes obvious. However, in the non- 
linear sigma models, one uses another regulator. Namely, since we have taken the 
inner product between x-eigenstates as (x\y) = g(x)~ l l 2 5(x — y), the momentum 
operator is represented by pi = —hig'^^^dig 1 ^, and this representation is clearly 
obtained from pi = —hidi by a similarity transformation with g~ 1//4 (x). As a result, 
the regulator used in non-linear sigma models is no longer g^^d^^fgg^dyg^^ but 
rather g^^d^^g^dy. In terms of momenta this reads g~ 1 ^Pig 1 ^ 2 g l -'Pjg~ 1 ^ 4 ■ We 
now explain in more detail how the Einstein anomaly vanishes if we take <fi as an 
independent field. 

Given a set of symmetries one wants to preserve at the quantum level in field the- 
ories, there exists a method to construct a regulator which yields consistent anomalies 
and which preserves these symmetries |18| . The basic idea is to construct a mass term 
which is bilinear in fields and which separately preserves the symmetries. For scalars, 
an Einstein invariant mass term is clearly M00 with = g 1//4 0. The regulator is then 
the kinetic operator for these fields, i.e., R = g~ l ^ A d^ L g l / 2 g IJ,u d u g^ 1 ^ . The Jacobian 
for Einstein transformations reads 

= i + ed, + \{o,e) = i + + d„e) (m 
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and using an orthonormal basis 4>n(%) satisfying / <j)* M (x)(j)N(x)dx = 5mn the anomaly 
becomes 

An = 4>* N (x)(eo ll + d ll e)Mx)dx 

= \\ & N (x)ed^ N (x)dx - ~J d^* N (x)CMx)dx. (104) 

It is obvious that this vanishes if the basis {4>*n} is in one-to-one correspondence with 
the basis {4>n}, as is the case for e.g. plane waves. In a more invariant language, 
the vanishing of the Einstein anomaly follows from the fact that the Jacobian is real 
and anti-hermitian, and the regulator is real and hermitian. The trace of the product 
of two such operators always vanishes, which can be deduced from the fact that the 
trace of any real operator (i.e. an operator that satisfies A{<p*) = A(<f>)* ) is equal to 
the trace of its hermitian conjugate. 

However, as explained above, in order to evaluate anomalies by using equivalent 
non-linear sigma models, the basis is not directly compatible with the conventions 
for (x\y) we have chosen. For chiral anomalies or trace anomalies, it does not matter 
whether one uses <fi or 0, because a similarity transformation with g 1//4 has no effect 
on 75 or 1: g 1//4 75g~ 1//4 = 75 and g 1 l i lg~ 1 l A = 1. But for gravitational anomalies, 
the basis <ft is more convenient. As follows from appendix A. 3, one finds for Einstein 
anomalies (= gravitational anomalies for <fi) in a covariant notation 

J = eD,+ l -(D,e) (105) 

where = — |r M! / and (D^) = + ry„£". The Y ^ terms cancel in J. 
The regulator is 

R = g-^d^g^dvg- 1 / 4 = g^R^g-V* (106) 

with R cov = g~ 1 ^ 2 d^ j g 1 ^ 2 g^ u d u the usual covariant scalar D'alembertian. Using cyclic- 
ity of the trace, or equivalently, making a similarity transformation (change of basis), 
the Einstein anomaly becomes 



A n (Em) = Tr (jexp(.R/M s 



= ^Tr (g- l '\eD» + D^g 1 '* exp(R cov /M 2 )) (107) 

A direct evaluation of this trace using plane waves is given in [37]. In the non-linear 
sigma model, jd^ becomes g 1 ^Pig~ l ^ A , and thus 

An = ^Tr(?pi + fig) exp(-(3H/h)) (108) 
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with H = g~ 1 ^ 4 Pig 1 / 2 g l: 'pjg~ 1 / 41 . Note that J comes out Weyl ordered, so that ^4 n (Ein) 
can be evaluated using the methods of the preceding sections. This explains why H 
is the regulator. (Of course, the Einstein anomaly vanishes, as we already explained, 
but other anomalies can be calculated with the same regulator.) 

For spin- 1/2 fields, the choice x a = x a g 1 ^ 4 as basic field leads to the same Ein- 
stein Jacobian, the field operator is now g x ^Jj) g^ 1 ^ 4 (where Ip is the usual covariant 
Dirac operator), and, as explained in ||18|| , the consistent regulator which leads to 
vanishing Einstein anomalies is now the square, namely R = g l ^Ip ij) g' 1 ^. Gravi- 
tational anomalies are really local Lorentz anomalies, but by taking a suitable linear 
combination of Einstein and local Lorentz transformations one obtains covariantly 
looking expressions. For chiral spin- 1/2 fields, one obtains then (see appendix A. 3) 

X(grav, spin-1/2) = -1tt(^ + D£») exp(i?/M 2 ) (109) 

where _D M = <9 M + ^cu At ab 7a7£> — \^^ v \ but the Y ^ terms in the Jacobian cancel again. 
Making the similarity transformation with g 1 ^ one finds (see ( |111|) ) 

A(grav, spin-1/2) = -^Tr^V, + ir£) exp{-pH/h) (110) 

where 7Tj = pi — ^ihui a bip a "4 !b and H given in (|52|) . This explains why we took this 
particular quantum Hamiltonian in the non-linear sigma models. 

Finally, for spin-3/2 similar results hold. However, here the Jacobian is more 
complicated. In an expression is given which does not correspond to a linear 
combination of local symmetries of supergravity. However, as we explain in appendix 
A. 3, the difference (which would be difficult to evaluate for non-linear sigma models) 
vanishes if one uses as regulator the same regulator Ip JJ) as for spin-1/2. Why this 
is the correct regulator is also explained in appendix A. 3. 

3.1 Chiral anomalies in gravitational couplings 

The simplest anomaly one can calculate by means of the path integral methods we 
have developed, is the 75 anomaly due to a loop of a spin-1/2 field coupled to ex- 
ternal gravitational fields. As regulator for spin-1/2 fermions we take JJ) Ip , which 
can be rewritten as the sum of a d'Alembertian and a gravitational curvature term 
(Weitzenbock identity) 

Iplp = ^A + ^VtA,^] 
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= ^D^^Df + R l3ab ^)l a l b 

= ±D?>Jgg*Df> + \R (111) 

where D'A = g ij DiDj = g i i(Df ) D j - T^D k ) and £>{ 0) = d t + \uj ia bTl b '■ We used 
the cyclic identity Ri\j a b] = to replace ^^ a ^ b by e a ^ b — e b ^ a , and then we used 
that the Ricci tensor R ia is symmetric. 

In the non-linear sigma model we choose the representation 

di - ^ V W 1/4 ; 7*e? = 7 a - (112) 
Then {ip a ,^ b } = 5 a6 , and the Hamiltonian becomes 

- yoWV 6 ) W (ft " f ^ - ^ (113) 

The representation p^ = g~ 1 ^j-^jg 1 ^ is fixed by the requirement that pi be hermitian 

and that the inner product is given by (x\y) = —?= — S(x — y), from which we find the 

v»( x ) 




completeness relations Jd n x Jg(x)\x)(x\ = 1 = Jd n p\p)(p\, and the inner product 
(x\p) = {2'Kh)~ n / 2 exp • xj g^ 1 / 4 . In the Fujikawa approach to anomalies, one 
often uses plane waves to evaluate traces in field theory. Since these plane waves 



are normalized to / exp(jrp(x — y))d n p = (2irfi) n ' 2 5(x — y), the regulator is in these 
cases exp(— Ip Jj) /M 2 ) with Ip = g 1 l A lfi g^ 1 ^- In the non-linear sigma model we have 



inner products (x\y) = (l/yg(x))S(x — y), and now the regulator corresponds to 
exp(- J J /M 2 ). 

We recall that under Einstein transformations with anti-hermitian generator E = 
^{pj,F j (x)} for the orbital part, x % — »■ x l + [x\E\ = x l + F l (x) = x'\ while the 
momenta transform as pi — > p\ = \ \^^iiPj\- Clearly ip a does not transform if we 
take x\ pi and i\) a as independent variables. It follows that also m — Pi — ^-uj iab il) a il) b 
transforms as t\\ = llf^r? 7 ^} (after adding a spin part to E which completes the 
transformation of Ui ab to that of a vector), and this shows, in the same way as in 
the bosonic case, that H is Einstein invariant. Since pi is hermitian, H is clearly 
hermitian, too. 

In a similar manner we may construct the orbital part of the Lorentz generator 
L = ■^■\ bc (x)^ b ^ c and show that 8lPj = ■^(dj\ bc )^ b ip c , 5i^> a = ^\ a c *p c and so, after 
adding a spin part to L, 5 L (uj iab il) a il) b ) = —(d i X ab )'ijj a 'ilj b , from which also the Lorentz 
invariance of H follows. 
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In the non-linear sigma model, H = ^QQ, where Q = \/2^ a e l a g 1 ^'Kig~ 1 ^ is the 
super symmetry generator. Since [H, Q] = 0, this Hamiltonian is supersymmetric. 
One can easily verify that Q is Einstein Lorentz invariant. 

The chiral anomaly is given by 

A n = Tr 75 e-^. (114) 

We shall compute this expression both in the case we double the number of Majorana 
fermions and in the case we combine the Majorana fermions into half as many Dirac 
fermions. 

We start by doubling the number of fermions, in which case we evaluate the trace 
in an artificially extended Hilbert space. After we introduce the free fermions ip% 
(a = l...n), we have (for even n) 2 n / 2 states in the ipi sector and 2 n / 2 states in 
the ip2 sector (combining the n ipl and n ip% in n pairs of creation and absorption 
operators). Hence, we must divide the trace over ipi and tp2 by a factor 2 ra//2 , since we 
really should only take the trace in the ipi sector. 

We shall now first express 75 into ip a . We define 75 by 

75 = {-i) n/2 lil2 • • • 7n n even (115) 

So, in n = 2 one has 75 = r 3 , and in n = 4 one has 75 = —71727374- In general 7 | = 1 
and 75 is hermitian. Identifying 7" = a/2-^i = + a — 1 ... n, we can evaluate 
the matrix element of 75 between fermionic coherent states and find 

n n 

(z\k\v) = H) n/2 m ii(v a +n = H) n/2 U(v a +n (m 

a=l a=l 

The expression for the anomaly becomes (recall the factor of 2 n / 2 ) 
An = 2~ n ' 2 [ f[dxiJgJx~oj [ ft dfj a drj a dCdC 

J i=l J a=l 

e^lslv) e~ m (x , v\e-^\x , (117) 

The last term in the above expression, the transition element, contains a factor e^, 
and further contributions from loops (which depend only on the sum (£ a + f] a )), but 
no classical action since the trace puts the initial and final points in x-space equal to 
each other. Now write 

n n 
a=l a=l 
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and perform the integral over 77 and £ (rewrite the measure d^ a dr] a in terms of the 
variables 77 — £ and 77 + £ as 2 n <i(£ a + r] a )d(r] a — £ a )). We find 



•An 



(_i)n/2 



i=l a=l a=l 

n n 

iid4vw) / n#b g e" iSioop8( "° , ^ ) 



(2vr/%)™/ 2 7 ^^V^^j ^r^^ (U9) 

We were able to integrate out the -02) canceling the factor 2~ n//2 , because Si oops only 
depends on the combination ip£ g = -^j(£ a + ?7 a ), which is the background value of -0i- 
The factor (2irhf3) n / 2 normalizes the leading (classical) singularity in (x\ exp(— |if)|x) 
to a Dirac delta function. The loop contributions Si oops (x , ip^ g ) are defined by 

M~^s lO o PS (x ,^ e ) = / P -^s int \ ^20) 



with the propagators given in fl80|). S- mt is the interaction part of the action S 

q 

S = -M a ^ a (0)/(0) + [dt 








/* 




-0 





+ M ab ip a ilj + -x'tUiabi/j^l + (ww and IT terms) 



121^ 



with the fields subject to the boundary conditions x l (0) = x l {— 1) = Xq and ^"(0) = 
fj a ,ip a (—l) = £ a . We now rescale t = /3t, which leads in the bosonic sector to 

1 

(122) 



h bos (3h 



dr -gij(x) (x l x j + b l c j + a 1 a? 



To obtain also a factor -| in front of the fermionic terms we rescale ip and ip suitably 



Then 



1 9 _1 



1 , 



dr [ 5 ab ^' b + -xW^'V 



1 



(123) 



W a (0)^ ,b (0) (124) 



The same rescaling is applied to the corresponding background values 



X" 



iX 



fa 



vp' 1 ' A (125) 

With this rescaling of the fermionic background fields, the (3 dependence in the mea- 
sure in ( |119[ ) is also canceled. Dropping the primes, we arrive at 



A 



i=l J a=l 



(126) 
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with the action the sum of ( 122 ) and ( 124 ). The IT and lulu term are of order (3 (due 



to the rescaling t = f3r) and play no further role in the evaluation of chiral anomalies. 
Expanding 

x* = 4 + q\ r = X a + < u > r = V a + < u (127) 

we see that 

(i) all vertices are proportional to 4 and all propagators are proportional to (3. 

(ii) hence only one- loop graphs contribute. 

(hi) In a frame where Ldi a b{%o) — 0, there are no terms linear in the quantum fields 
in S, so no tadpoles. Therefore, only vertices with exactly two quantum fields 
are relevant for the one-loop graphs in that case. 

(iv) Consequently, by expanding gij(x) and uii a b{x) around xq, and after substituting 
( |127f) in S, we find that in the frame Ui a b(xo) = the only vertices are 

1 1 / 

S^t = ^{djCOiabixo)) / dr gVCA (128) 



which can be rewritten as 



-l 



1 1 



h S int = —Ri^uix^^Jdrq^ (129) 

Thus the one-loop contributions are due to a q loop, with at each vertex two 
sticking out. 

The one-loop result can now easily be evaluated by noting that it is proportional 
to the one-loop determinant 

-1/2 



(_ j)n/2 r n 

An = J dUx ° v 9{xo) n d w* 



a=l 



det [—^9ij(x ) + iRijabKg^- 



b _d_ 

dr 



det (-^flfy(xo) 



(130) 

The denominator in this expression is the ghost contribution to A n . The factor % n l 2 
can be removed by a rescaling of ip^ g , in which case the only thing that changes is that 
Rijab is replaced by %Rij a b- With this rescaling the operator in the numerator of (|130|) 



is hermitian and it becomes manifest that A n is real. From here on, one can follow 
Alvarez-Gaume and Witten 0. First one skew diagonalizes Rij, then one computes 
the eigenvalues of the operator in the determinant. The product of the corresponding 
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eigenvalues then immediately yields the well-known result for the chiral anomaly as 
the A-genus of the manifold. An alternative derivation can be found in appendix 
A. 4, where the one-loop diagrams contributing to the one-loop determinant ( |130| ) are 
explicitly evaluated. 

We will now evaluate the trace in ( |114j ) in case we do not double the number 
of fermions, but instead combine the Majorana fermions into half as many Dirac 
fermions. In this case the expression corresponding to 75 becomes 

n/2 1 n/2 

75 = (_i)«/2 TJ + q A )-^ A - $ A ) = JJ (1 - 2f A ^ A ) (131) 

A=l 1 A=l 

where the ^ A and ^> A are defined as in (|54|) and we identify again 7° = \/2ip a . Since 
Pa = ty A ^ A for fixed A is a projection operator (P| = Pa) we can rewrite this as 

n/2 n/2 

7 5 = [] e" i7rVl/ * = exp ( -in £ ^ A ^ A ) (132) 

since expansion gives IIa(1 + (e~ i7r — 1)^/ A ^ A ). 

The matrix element of 75 between coherent states simplifies even further 

n/2 n/2 

(flTsto) = & II (i - 2 ^V) = n a - 2eV)(i + eV) 

A=l A=l 

n/2 _ n/2 _ 

= n(WV)=exp -£eV) (133) 

A=l V A=l ' 

To evaluate the anomaly, we put together the expression for the Jacobian and 
the transition element 

n/2 



A n = /n^o^S) / I[df) A dri A dt A dt A 

J i=l J A=l 

e«(e|75|?7}e-^(xo,77|e-^|xo,0 (134) 
As before, we extract a factor exp?y£ from the transition element, and write 

in „ n/2 

An = (27r/m) n/ 2 / UfcoyfiiM J T[ dfj A dr) A d£ A d£ A 

e tt-iv-vv+m exp (-ls loops (x , V A , e A )) (135) 

We can now trivially do the integral over r\ A and £ A , after which the above expression 
equals (|119[ ), with the identifications 

r hs = -L(^)/ 2 + ^)/ 2 ) a odd 

K g = -^(C /2 -V a/2 ) a even (136) 
Hence, also in this case, we obtain the same expression for the chiral anomaly. 
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3.2 Chiral anomalies in Yang-Mills couplings 



We consider complex spin-1/2 fermions coupled to external Yang-Mills gauge fields 
corresponding to a group G, with the fermions transforming under a representation 
R. We will only consider flat space here; by combining the techniques in this and 
the previous section, one can easily obtain the combined gravitational and Yang- 
Mills anomalies. We leave this as an exercise to the reader. New in this section is 
the treatment of the Yang-Mills ghosts by path integrals. In 0, they were kept as 
operators. 

The Jacobian is still 75, but the regulator is now proportional to 

ipip = d^d, + iyviA*: D A, ( 137 ) 

where = <9 M - gA*T a , [T a , Tp] = f a pT^ and [Dp, D v \ = -gF* u T a . Hence, com- 
pared to the bosonic case, there is now an extra term proportional to the Yang-Mills 
curvature. We represent the matrices (T a ) M N , M, N = 1 . . . dimR, where dimR is the 
dimension of the representation R, by operators 

f a = c* M (T a ) M N c N , (138) 

and require the anticommuting relations 

{c N ,c* M } = 6 N M . (139) 

Then [T a ,Tp] = f a pT 7 . The Dirac matrices 7° are again represented by y/2ip a with 

{iP a ,^ b } = 5 ab (a,b=l,...,n) (140) 

We will only represent the Majorana fermions ip a by half as many Dirac fermions 
here; if one doubles the number of fermions one obtains the same answer. This can 
easily be demonstrated in a similar fashion as we did in the previous section. With 
half as many Dirac fermions we had ( |132j ) 

7 5 -> e ~ l7r ^ A . (141) 

The anomaly is then represented in the non-linear sigma model by 

An = TrV^^V^ (142) 

1 h 2 

H = -( Pj + MA^c*T a c)(p k + mA^c*Tpc)^ k + -r^ b F: b c*T a c (143) 
I J z 

We have rescaled A"(x) and such that the coupling constants have been absorbed. 

40 



The prime on the trace, Tr', indicates that the trace is not over all states in 
the Fock space (which are |0), cfojO), <? Ml (? Ma \0), ■■■,c* Ml --- c* MdiinR \0), but rather only 
over the one-particle states c* M |0). Only on these states does c*T a c act like the matrix 
T a . To still write the trace as an unconstrained trace, we introduce the one-particle 
projection operator P. We claim that 

P =: xe~ x :, x = c* M c M . (144) 

Indeed, from the definition of x it follows that 

: x n := (x-n + 1) : x n ~ x :=••- = n\ ^ X j . (145) 

Then we use a representation of the Kronecker delta S X) i, 

d 

5 X i = limxp(l - p) x ~ l = limp— [-(1 - p) x ]. (146) 
Expanding (1 — p) x in a power series and performing the limit p — > 1 we arrive at 

EnC-l^f X ) = E^f (147) 

n>0 V " / n>0 V " / ( ^ 

Using ( |145| ) we get the desired result, namely 

P = : xe -» : = 5 X>1 . (148) 

The anomaly is thus given by 

A n = Tr (e~^ A : c^e"^ : e^ H ' h ) . (149) 

Using complete sets of coherent states, one finds the corresponding path integral 
representation 

A n = tr x tT f tr gh (xo,x g h,Xf\e- l ^PI gh Ife- mh \xf,X 9 h,x ), (150) 

where 

Igh = J df] gh dri gh \ri gh )e~ n ^ h (f] gh l (151) 

I f = jdf} f drj f \rj f )e-^f(f} f \, (152) 

tr x . = J dx (153) 

tegh = J dx g hdx gh e* shXsh , (154) 

tr/ = Jd X fdx f e*' x '. (155) 
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We have used that in fermionic spaces the trace of an operator is given by ( p2|) and 
further inserted two unit operators / for which we used the decomposition in terms 
of coherent states. 

The trace involving e~ l7T ^P factorizes. The ghost dependent part of A n reads 
dXghdx 9 hdf] gh dr] gh e^ hX ^e~^ h ^ h {xgh\ ■ xe~ x : \Vgh)(v g h\e~ PH \Xgh)- (156) 



Since P is a one-particle projection operator its matrix element simply yields 

(Xgh\ ■ xe~ x : \r} gh ) = XghVgh, (157) 
The integration over x gh then yields 

/c^i(EwS^ = z(( n xi)< h { n xi)) 

J M M N>M N<M 

= Exf R -4-i (158) 



M 

namely, a product of all x^h except that the M-th factor is replaced by r^. The 
integration over rj gh then yields 

/ n n *&)) =e n (159) 

J M N>M N<M M N^M 

Clearly, this operator projects an arbitrary function of fj g h and Xgh onto the terms 
with precisely one f] gh and one Xgh- 

In the fermionic sector, the matrix element of 75 is again (|133|) 

(Xf\e- in ^ A \Vf) = e-* fVf . (160) 



This leads to the integral 

(l\f(lijf( x ' x '< *ns t = e -f? W. (161) 



Inserting these subresults from the ghost and fermionic sectors into the anomaly 
equation we arrive at 



An = J II [dXo d X%h d Vgh,M d Vf,Adx J j 



eft 

i,M,A 

X [E( II Vgh^Xg^y^ixo^^Vjle-^lxf^gh^o). (162) 

M N^M 



The transition element contains a factor e' ?/X - f and a factor e V9hX3h , in addition to 
contributions from loops. The former cancels against ( |161|) , so that we obtain 

d,X l n 



oVghXgh 



j y/^irnp j M ^ A M N _^ M 

x exp ^-^Si oops (x , fjg h , fj f , Xgh, X/)) , (163) 
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where 

e nfXf+nghXgh-\Si aaps {x ,n gh ,fjf,X gh ,Xf) _ /g-sSint^ (164) 

Here, S- m t is the interaction part of the action 

° dttM^WV - ^r^FZtclATc)" „c N ), (165) 

The couplings xAc*Tc result from integrating out the momenta p, and combine with 
the ghost kinetic term to the covariant derivative D t c M = c M — x J A^{x)T a . The 
fields satisfy the boundary conditions x(0) = x(—/3) = x , $/(0) = fjf, c*(0) = fj gh , 
if)f(—/3) = Xfi an d c(—(3) = Xgh- Decomposing all fields into background parts and 
quantum fluctuations as in ( p.27|) , we note that the object Si 00 ps now contains tree 
graphs because the four-fermion couplings ipx c * c contain terms linear in quantum 
fields and hence lead to tree graphs. These tree graphs do, of course, contribute to 
the classical action, as we discussed in section 2. 

In order to study the loop expansion in detail, we rescale t = /3t, and if) — > 
ip'iph)- 1 / 2 , if; -> V'(^)" 1/2 , Xf -> Xf'W' 112 , fjf -> fj/iph)- 1 / 2 . After the rescaling 
of dxf and dfjf the measure becomes completely f3h-independent, and 

2 



+ £ dT(c;,«i M - i^cMrj'V + iv>'-FSc^(r a ) M K c iv ) (166) 



Since the a; and ?/> propagators are proportional to (3U and all vertices and c*c 
propagators are (3Ti- independent, we need not consider the vertices containing the 
<?\ ^qu,A or ipq n of (|127|) . Hence, we can restrict our attention to the vertices Fipipc*Tc, 
with the ^'s replaced by their background value. There are no classical contributions 
to S^ops since we evaluate the classical action from y to z with y = z = xo- So 
altogether we are left with 

II [^= d X g tdf] g h,M d Vf,A d xf}[Y,( II %h,NXgh 
i,M,A V 27171 M NjtM 

1 1\ 

(167) 



oVghXgh 



exp 



^ dr(c^c M + -^ g <F a a 6 c*r Q c) 



loops 



where V'bg are the classical values of ?/> a exactly as in ( |136| ). 

The propagator (c M (cr)c* N (r)) is equal to 5^9(a — r), hence no closed c-loops 
can contribute (in a closed c-loop one always moves somewhere backwards in time). 
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Only terms with precisely one fj g h and one Xgh can contribute due to the projection 

operator in the measure. These terms are just tree graphs, with c* = fj g h at one end 

and c = Xgh at the other end. If we consider a diagram with k vertices, then it will 

contain an integration 

r° 1 
J ddf- da k 9(a l - a 2 ) ■ • ■ Q{o k -i ~ = jj- (168) 

The combinatorical factor l/k\ from the expansion of e~ s l n cancels because there are 
exactly k\ different ways to build a connected tree graph from the vertices. Putting 
all this together yields 

- ^icops = fj gh>M (e- - if NX » h (169) 
M 



The background term exp(f] gh:M x g h) in ( |167|) cancels the —1 in ( |169|) . The integration 
over fj gh and Xgh in (|167|) with «Si oops replaced by ( |169| ) can now easily been done, since 
it picks out the piece with only one fj g h and Xgh from S^ops- Transforming in addition 
variables from f]f and Xf to leaves us then finally with the following result 

• \ n/2 „ 



n/2 1 



J dxie a ^tr(F aia2 ---F an _ ian ) (170) 



AnJ (f)! 

where the trace is over the Yang-Mills indices of T a in = F^ b T a . This is the correct 
anomaly. 



3.3 Trace anomalies 

We can now easily compute the trace anomalies for a spin-0 and spin-| field in an n- 
dimensional quantum field theory. In this case, the Jacobian is equal to one, so we just 
need to evaluate the trace of the appropriate transition element, and no additional 
operators need to be inserted. The anomaly is then equal to the /9-independent part 
of the trace. New in this section is the treatment of fermions by path integrals; in 
they were kept as operators. 

For a real scalar field in n = 2 dimensions, we can directly take the trace in flTT|). 
Singling out the /3-independent part yields 

^trace = ^ R (171) 

2 24tt v ; 

In order to obtain the trace anomaly for a scalar field in higher dimensions, we need to 
compute the terms of higher order in /3 in the transition element. Since the transition 
element contains a factor (27r/3/i)~ n//2 , one needs f + 1 loops in n-dimensional space. 
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We will now evaluate the trace anomaly for a spin-^ field in n dimensions. In this 
case we have to compute the trace of the N=l supersymmetric transition element, and 
project out the /5-independent part. As claimed before, both approaches to obtain 
fermionic creation and annihilation operators will lead to the same result. 

Recall that when we double the number of fermions we should normalize the 
trace by dividing by an extra factor 2 n//2 . The trace anomaly for a spin-| field in n 
dimensions is therefore given by (cf. 0, equation (2.9)) 

Xn m 1 = -g^2 timJ d X a dfj a e r]X (x ,fj\exp (~^J \ x o,x) ! a = l...n 

(172) 

with the transition element given in (|8l|). For n = 2 the trace anomaly becomes 
which is indeed the result for a Dirac fermion; for the anomaly of a Majorana fermion 
we have to divide this expression by two. 

When we instead combine the Majorana fermions into half as many Dirac fermions, 
we should take the trace of the transition element given in (|96|). In this case of course 
no extra normalization factor is needed, and we find directly 

AT n ~^ = -hlimJdx A df} A e f > x (x ,f)\exv \ x o,x) ! A = l...n/2 

(173) 

now with the transition element given in fl96|). This gives again the same result for 
the anomaly as above. 



4 Conclusions 

In this article we have given a complete, explicit derivation of quantum mechani- 
cal path integrals for bosons and fermions, both for Dirac fermions and Majorana 
fermions. Our main result is that the factors 5(a — r) in the Feynman rules for con- 
figuration space path integrals should be interpreted as Kronecker delta functions, 
even in the continuum case, and should not be regulated by mode regularization. 

We define our path integrals in curved space by starting from the Hamiltonian 
(operator) formalism. After inserting complete sets of states (coherent states for 
fermions), and Weyl-ordering the Hamiltonian (leading to order h and ft 2 terms in the 
path integral action), we obtained the discretized propagators and vertices in closed 
form. With these ingredients one can construct a loop expansion of the path integral 
in terms of Feyman integrals. Some of the bosonic propagators contain 5 (a — r) 
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singularities, but adding 'Lee- Yang ghosts', terms with two or more 8{a — r) cancel. 
Terms with one 5 (a — r) should then be evaluated as indicated above. 

We paid particular attention to Majorana fermions. Of course, starting with an 
arbitrary initial state \A) and acting on it with products of "0's, the states so obtained 
will span a Hilbert space on which the r/>'s can be represented as matrices. One 
can then define matrix- valued Hamiltonians as in ||. We found it much simpler to 
define creation and annihilation operators and then to use the standard formulation 
of coherent states pifl. The particular way of defining creation and annihilation 



operators is, of course, arbitrary, and so is therefore the choice of vacuum, but in 
problems involving a trace over the Hilbert space, this arbitrariness should cancel. 
We achieved the construction of creation and annihilation operators in two ways: 
either by combining the Majorana spinors pairwise into creation and annihilation 
operators ('halving'), or by adding another set of Majorana spinors ('doubling', this 
works also for odd-dimensional spaces). Of course, the Hilbert spaces, vacua etc. are 
different in both cases, and indeed we found different expressions for the transition 
element, but the anomalies came out the same. This confirms our claim that in traces 
over the Hilbert spaces, differences created by choosing different vacua should cancel. 
We verified our formalism by computing the transition elements for a bosonic and 
several fermionic transition elements through order /3, and comparing the results with 
those from (unambiguous) operator calculations. We provided further evidence by 
doing a three- loop calculation in appendix A.l. 

We applied our general formalism to trace anomalies and to chiral anomalies for 
spin-1/2 fields. The chiral anomalies were already studied by Alvarez-Gaume and 
Witten and the trace anomalies in ||, but we have treated the Yang-Mills ghosts 
and the Majorana fermions on equal footing with the bosons and obtained a uniform 
path integral treatment. Moreover, our derivation makes a detailed and complete 
treatment for any anomaly possible, including all normalizations. 

On the more technical side, we have seen that the Hamiltonian H in the operator 
formalism and the Hamiltonian function H in the path integral are related by the 
formula H = (H)w, where (H)w contains in general terms of order % and Ti 2 . If H 
is Einstein invariant, the corresponding action in the path integral is not Einstein in- 
variant (although the transition element is); rather there are noncovariant IT terms. 
Conversely, if the action is the naive action, H will contain K and h 2 terms. In par- 
ticular, the supersymmetric Hamiltonians (whose ambiguity was fixed by requiring 
hermiticity and Einstein invariance) do not lead to the usual classically supersym- 
metric action in the path integral. Rather, there are extra terms proportional to h 2 . 
For chiral anomalies, though, these extra terms in the action do not contribute, which 
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explains why the results of Alvarez- Gaume and Witten are correct (that they are 
correct can be checked by doing loop calculations in the corresponding quantum field 
theory, see and [@). For trace anomalies, the extra terms do matter. Here higher 
loop calculations are needed and we stressed that the noncovariant vertices (of the 
form IT and ujuo) as well as our new Feynman rules must be taken into account, even 
when one uses normal co-ordinates, to obtain the correct results. 

This concludes our analysis of quantum mechanical path integrals. One might 
wonder whether the subtleties we have found in one dimension have a counterpart in 
higher dimensions. For local field theories this seems unlikely, because we expect that 
possible extra terms in the path integral will be proportional to 5 n_1 (x), and hence 
would vanish in dimensional regularization. However, in non-local field theories, such 



as Yang-Mills theories in the Coulomb gauge ]T3|, |I5 , there might be effects. This is 
under study. 
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A Appendices 

A.l A 3-loop computation 

We check that there are no counterterms beyond the two-loop counterterm in (|T3J) by 
performing a 3-loop calculation in a model which describes a free particle, but which 
has been cast into the form of a non-linear sigma model by a nontrivial co-ordinate 
transformation. This extends the two-loop phase space calculation of |LJJ. We shall 



then redo the calculation in configuration space and show that we obtain the same 
result (Matthews' theorem). 

We consider a free massive point particle on the interval — oo < t < oo with 
L = ig 2 — \q 2 and substitute q = Q + |Q 3 - Then the action becomes 

HQ, Q) = \q\i + Q 2 ) 2 - \Q 2 {i + \Q 2 ) 2 (174) 

and the Hamiltonian is H(Q, P) = \P 2 (l + Q 2 )~ 2 + \Q 2 (l + \Q 2 ) 2 . The first-order 
action in phase space reads 

L(Q,P) = (PQ-]-P 2 -±Q 2 ) + £ mt (Q,P) 
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L in t(Q,p) = ^ 2 (2Q 2 + g 4 )(i + g 2 )- 2 -ig 4 -^g 6 

= p 2 Q 2 -\p 2 Q a -\q a -^Q & + ... (175) 

We consider the vacuum self-energy at the 3-loop level. There are 

(i) "clover-leaf graphs" from the P 2 Q A and Q 6 couplings with three loops all meet- 
ing at one point, 

(ii) a string of three loops with two four-point vertices. We call such a string a PP 
string (or PQ string or QQ string) if the equal-time contractions at the ends 
of the string consist of two P lines (or one P and one Q line, or two Q lines. 
Since an equal-time PQ contraction vanishes, this exhausts all possibilities.) 

(hi) watermelon graphs with 4 propagators between 2 vertices. 

(iv) a one-loop contribution from the extra two-loop interaction AV; since AV 
is of order h 2 , this one-loop graph contributes also at order h 3 . To evaluate 
= \T) k T j u g kl we note that = (1 + Q 2 ) 2 and find 

AV = ^Q 2 (1 + Q 2 )-* = ^Q 2 + ... (176) 

There is no R term since the metric is flat in this model. For the phase-space 
calculation we use the propagators 

<P(a)P(r)> = <Q(a)Q(r)>= 1 -e-^ 

<P(a)Q(r)> = - < Q(a)P(r) >= ZV'I— l e (o- - T ) (177) 

and further the equal-time contractions 

< P{cr)P{a) >=< Q{a)Q(a) >= ^; < P(a)Q(a) >=< Q(a)P(a) >= 

(178) 

48 



We find the following results 



From P 2 Q 2 - %P 2 Q 4 



clover = — 



9i 
16 



From l „Q A l*Q G Cross terms 



clover = — 



5i 
18 



16 



P — P string 
Q - Q string = ^ 
P — Q string = 
watermelon (no PQ) = ^ 
watermelon (two PQ) = \ 
watermelon (four PQ) = ^ 



string 



watermelon = ^ 



QP string = -f 
QQ string = | 

watermelon (two PQ) 



From AV : 



(179) 

Adding all contributions, we find the correct result: the sum of the vacuum self- 
energies vanishes at the three-loop level. This demonstrates that at the 3-loop level 
no further counterterms are present in the phase space path integral. 

In the configuration space path integral, there are extra vertices and an extra term 
in the QQ propagator. According to Matthews' theorem, the final answer should be 
the same as in the phase space approach. We now check this. The QQ propagator 
is the same, while the QQ propagator in configuration space is equal to the PQ 
propagator in phase space. The QQ propagator differs from the PP propagator by a 
Dirac delta function 



l 



-i\a— t\ 



+ i5(a - r), 



;iso) 



which is due to differentiation of the time-ordering 6(<j — t) and 8(t — a) functions. 
The action contains also ghosts, 

HQ, Q, b, c, a) = ^(QQ + bc + oo)(l + Q 2 ) 2 - \q 2 (1 + ^Q 2 ) 2 , (181) 

whose only role is to cancel products of Dirac delta functions. We shall now show 
that the contribution from the extra vertices in the Hamiltonian approach equals the 
contributions of the extra terms with one Dirac delta function in the configuration 
space approach. 



Q=P 



The extra vertices in the phase space approach are L int (P, Q) — L int (Q, Q) 
-2P 2 Q 4 . Hence we get only an extra contribution to the clover leaf graph. Since 
-3/2P 2 Q 4 gave — 9z/16 according to fll79D , we now get 4/3 times this contribution 



contribution extra vertices m phase space = — — . 
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;i82) 



In the configuration space approach, all delta functions in the clover leaf graph 
cancel, since for every QQ contraction there is a compensating be and aa contraction. 
In the graph with a string of three loops, no delta functions remain if one of the outer 
loops is a QQ contraction. The graph with one Q in an outer loop vanish. Hence 
only the graph in which the inner loop contains two QQ contractions yields an extra 
contribution. This contribution comes from the graph with two Q 2 Q 2 vertices. It 
contains a term proportional to 8{a — r) 2 which cancels against the contribution from 
similar graphs with an internal ghost loop. Only the terms with one delta function 
remain and one finds 

extra contribution to string = — — . (183) 

Finally there are the watermelon graphs. We must consider graphs with two QQ 
propagators and graphs with one QQ propagator. However, the latter do not con- 
tribute because they also contain a (Q(a)Q(r)) and a (Q(a)Q(r)) propagator, and the 
integral / 5 (a — r)e(a — r)e(r — o)dodr vanishes according to our rules. In the graph 
with two (QQ) propagators we again need the term proportional to one 5 (a — r). 

extra contribution watermelon graphs = — . (184) 
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Since (|182|) equals the sum of (|183 ) and ( |184j) , Matthews' theorem is verified in this 



model at the three-loop level. The extra potential in ( 176|) , being proportional to 



two Christoffel symbols, will not contribute below the three-loop level in normal co- 
ordinates. However, at the three-loop level it does contribute. Note that the usual 
co-ordinate transformation from arbitrary co-ordinates to normal co-ordinates will 
yield correct results in the one- and two-loop computation, but will yield an incorrect 
result at the three-loop level. In our model, this is very clear: in normal co-ordinates, 
the action ( |174| ) reverts to a free model, but the extra potential in ( |176| ) would yield 
a spurious contribution in the normal co-ordinates q. 



A. 2 Higher derivative theories 

To test our rules in a higher- derivative and higher dimensional model, we consider a 
massive real scalar field (p with rather singular interactions |39 



£= -$d lt <pd ft <p-$m 2 <p i -$\(dp<p& l <p) 2 (185) 



One novelty that appears in this example is the need to introduce 'Lee- Yang' ghosts 
which couple to derivatives of the scalar fields, see (|203|) . 
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The conjugate momentum is denned by 

^ = S^x) = ^ + ( V? ' ^ tp,t ^) = ^ + X ^ (yvf-V 2 , (186) 

and the Hamiltonian density is given by 

oo 
n=0 

It is straightforward to find the first few terms of H, 

W<°> = I[^ + (V^) 2 ] + ^mV 

W* 1 ) = -£ m ,(^V^7r) = i[7r 2 -(V^) 2 ] 2 

^ (2) = i^^V^Tr)]^^ 2 ^-^) 2 ] 2 (187) 

The terms linear in A yield the same vertices as in the Lagrangian approach, namely 
£int{<Pi V<£>, 7r), but with replaced by 7r. Extra with respect to the Lagrangian 
approach are the vertices in TC^ 2 \ In the interaction picture where one uses the free 
field equations and in-fields ipi n and n in , one replaces i\i n by ip in . 

On an infinite ^-interval the boundary conditions on the bra and ket vacuum are 
that they are the lowest energy states. For the Feynman propagator one then obtains 

< 0\T M) \0 >= / gp^j^ (188) 
For the ipir propagator one finds 

< 0|7V(z)7r(y)|0 >=< 0\T V (x) <p(y)\0 >= A < 0\T V (x) V (y)\0 > (189) 

because hitting 9(x°—y°) produces a y°) times a commutator to); <f(y, to] 
which vanishes. However, for the ixix propagator one obtains an extra term 

< 0|T7r(V)7r(y)|0 >=< 0|Ty?(V)y%)|0 >= 

= dx^df < °I T ^^I° > S(x -y )Hx),7r(y)] 

< 0\T V (x) V (y)\0 > -ih5 4 (x - y) (190) 



dx° dy° 
More generally we have 



<0\Td li <p(x)dv<p(y)\0> = — —<0\T<p(x)<p(y)\0> 

-ih5,°5„°5 A (x-y) (191) 
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In momentum space we find 

< 0\T S M*)dMy)\0 >= / (0 ( ^ti, ~ «W*V) «»<->. (192) 

The same phase space propagators are obtained from the path integral by inverting 
the kinetic matrix of the fields ip and ir. 

We shall now consider elastic tptp scattering. There are the extra vertices and the 
extra terms in the unequal-time propagator to use, but we shall also need equal-time 
contractions. They are given by 

<0| Wl)Wl )|0>^^ = /||^ (193) 

These equal-time contractions are the same in the Lagrangian approach as in the 
Hamiltonian approach. We shall not try to regulate these divergent expressions, but 
we shall only use that for k^k u equal to kikj one may replace kikj by |<%/?, while 
for k^k v equal to kik one obtains zero. We shall denote the corresponding divergent 
integrals by 

p r d 3 k P 2 _ r d 3 k k 2 

I(kikj) = ^ tJ I(k 2 ), 1(1*1*) = (194) 

(where k$ = oo 2 ). The equal-time contractions in the Hamiltonian approach are thus 
equal to those in the Lagrangian approach. Our aim is to show that the extra terms 
with I{k 2 ) and /(/cq) m the Hamiltonian approach (due to extra vertices in H mt ) 
cancel algebraically with similar terms in the Lagrangian approach (due to the term 
i5®5®5 A (x — y) in the propagator). 

Recalling the vertices 

n mt = \n^ + \ 2 n^ + ... 

= ^d^d^) 2 + \\ 2 ^ 2 {d^d^) 2 + ... (195) 

there are 3 graphs to be computed, shown in figure 1. 

In (I) we find the combination (< p fM (x)p u (y) > +i<5°5°(x — y))(< P P {x)p a (y) > 
+i5° p 5®5 A (x — y)) where we have written the Feynman propagators as a sum of phase 
space propagators plus contact terms. The two cross terms yield extra terms pro- 
portional to the equal-time contractions < p^{x)p v {x) >. Since the phase space 
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I 



II 



III 



Figure 1: The three graphs to be computed. Arrows indicate where the noncovariant 
propagator or the extra vertex is to be used. 

propagators are continuous at x — y, there are no subtleties in taking the equal-time 
limit. In (II) the equal-time contraction is the same in phase-space and in configura- 



in the loop is the same in phase space as in configuration space for the same reason, 
but the propagator connecting the two vertices has a noncovariant piece. The net 
effect of these noncovariant pieces in the propagator is to contract the propagator to 
a point. Hence, all 3 graphs become topologically of the form of figure (II), which 
makes the cancellation possible. We now study in more detail whether the extra 
contributions from (I) and (III) in the Lagrangian approach are equal to the extra 
contribution from (II) in the Hamiltonian approach. 

For (I) we need the contractions in 



As first contraction we can take < d v ip(x)d a ip(y) >, it comes with a statistical fac- 
tor 4 x 4 = 16. The second contraction can then still be done in 3 different ways 
(d u (pd a (p, d v ipdpip or d p (pd p (p) ; however, to avoid double-counting, we need to multiply 
with a factor 1/2. Hence 



Si = ^-[{d^f{x)<dM^)dMy)><d v V {x)d^{y)>{d pV f{y) 

+ 4(d^)\x) < dM*)dMy) >< dr<p(x)d p <p(y) > d^(y)d p v(y) 

+ 4dMx)dvV>(x) < d^(x)d p ip(y) >< d u V (x)dMy) > d p <p(y)d a <p(y)} 



tion space, but the vertex Ti(2) is extra in phase space. Finally, in (III) the propagator 




(196) 



-A 2 16 



(197) 



In each pair of contractions we make the substitution 



< d li (p(x)d v (p(y) >< d p ip(x)d a ip(y) >= 

< P^)Pu(y) + - y) >< P P {x)pM + i%S°j\x -y)> (198) 
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and retain only the two cross terms, replacing S A (x — y)p p (x)p a (y) by equal-time 
contractions for which we then use (|193|). This yields 



Si = ~ l -{d^fl{kl)-2i{d^fd a ^I{k k a ) 

- 2id^d p (p(p 2 I(k fM k p ) (199) 



In Su we find four different contractions, leading to 



S 



11 



4tjfd IM <p&'<pI(k li k v ) + 2(p\dLpyi{k u k v ) 
+ ^d^(d^) 2 I(k^) + (dipyi(k , ko)] (200) 

Finally, in (III) there are two different contractions, yielding 

S IH = -i(d(p) 2 (p 2 I(k tl k") + 2i(d(p) 2 (p 2 I(k 2 ) (201) 

We tabulate the results by decomposing each (d^ip) 2 into (V<^) 2 — 2 , and for nota- 
tional simplicity we write k 2 instead of I(k 2 ), and k^ instead of I(k^). We find then 
the following extra terms 





(V^fk 2 


(V^k 2 


(Vif) 2 if 2 k 2 


(Vp)V*o ^ 2 


^k 2 


(/): 





-1/2 


-2/3 


1 + 2 


-1/2-2-2 . . 
' _ 3 (202) 


(III) : 








-1 


3 1 


(II): 





-1/2 


-2/3 - 1 


1+4+1 1 


-2-1-4-1/2 



Since the sum of (I) and (III) is equal to (II), we conclude that Matthews' theorem 
is satisfied in this example. Note that this conclusion is not based on a particular 
regularization scheme; rather the cancellation is algebraic. 

Finally we must account for the terms with [<5 4 (x— y)} 2 in ( 197 ). In the Lagrangian 
approach, one integrates out the momenta by a saddle-point method. This means 
that one must evaluate ^(pq - H(p, q)) = -^H(p } q) = —j^(pq(p) - L(q, q(p))). 
Using ^L(q, q(p)) = p^ one finds 

d 2 dq 1 1 

— [pq - H(p, q)} = - — = = -^j- 

U P U P 94 Qq2 

Hence the a, b, c ghosts one needs to add are given by 



, . . , d 2 L d 2 L 
L(ghosts) = b c + a- 



(dq) 2 (dq) 2 
In our case this yields 

L(ghosts) = b (l + A {(V(f) 2 - 30 2 }) c + same with a. (203) 
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The 8 A {x — y) 2 terms from ( 197 ) are easily seen to be proportional to [(Vy?) 2 — Sip 2 ] 2 , 
and the ghost loops cancel these singular terms since the be and aa propagators are 
proportional to <5 4 (x — y). 

This concludes our discussion of this example and of Matthews' theorem, which 
states that the Hamiltonian and Lagrangian approaches to perturbative quantum 
field theory are equivalent. We used it to test our rules in a one-dimensional and 
a four- dimensional model. Crucial was the correct definition of equal-time contrac- 
tions. It is often stated that equal-time contractions are ill-defined, but if one wants 
to compute loop corrections in non-linear a models, one must deal with them. Equal- 
time contractions are, in fact, already needed in a much better known area: the Ward 
identity for the self-energy of charged massive scalars coupled to photons are only sat- 
isfied if one includes equal-time contractions. These seagull-graphs can be computed 
with dimensional regularization and are nonvanishing, but since they correspond to 
< q(x)q(x) >, they are the same in Hamiltonian and Lagrangian formalism and are 
the limit of x tending to y of < q(x)q(y) >. The equal-time contractions we consid- 
ered included < q(x)q(x) > and < q(x)q(x) >, and these are not the same as obtained 
from the limit y tending to x. 



A. 3 The covariant spin-3/2 Jacob ian 15 

In order to obtain a covariant expression for the spin-3/2 transformation rule un- 
der space-time transformations, and as a consequence a covariant expression for the 
corresponding Jacobian, we must take certain linear combinations of Einstein trans- 
formations and local Lorentz transformations. For spin-1/2 fields, this is easy and 
well-known, as we now show. Afterwards we consider the more complicated case of 
spin-3/2 fields. To wet the appetite of the reader for this problem, we first quote 
the final results, which were obtained by Fujikawa for spin-1/2 fl7|l , used by Alvarez- 
Gaume and Witten |2| and further studied in |4(| £yj 

spin 1/2: 8 Aw j> = X ^D^ + ^(D^W 

spin 3/2 :<WV>m = x^D^m + -(D^Wm + [(V) - P"Xm)]i 

D mX n = e m »(d,x n + uJ» n P X p )iX n = e» n x» 

D^ m = d^lpm + UJ^i'n + ^Iplq^m ~ T^^^m (204) 

15 These results were obtained in collaboration with R. Endo. 
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The fields ip and ip m are world-scalar densities of weight 1/2, namely ip = g 1 ^ 4 ^ and 
VVt — g 1 ^^Pm with ip m = em^. Thus the covariant derivatives contain a term T ftv v . 
Further the spin-3/2 field ip m , the so-called gravitino field, is a Lorentz vector-spinor. 
This explains the term uj^ m n ip n in its covariant derivative. 

The last two terms in the spin-3/2 transformation rule describe a local Lorentz 
transformation which acts on the vector indices with parameter D^Xn], whereas the 
first two terms contain a local Lorentz transformation which acts on the spinor indices 
with parameter x M cJ At m „. At first sight these rules do not seem to correspond to a 
linear combination of the usual transformations. We now proceed to demystify these 
expressions. We begin with the easier case of spin-1/2. 

For spin-1/2, we define Einstein, local Lorentz and "covariant" transformations 

by 

= X^ + ^d^ 

SUx") = S E ( X n + Se L (xV n ) (205) 
Thus for the spinor field ip one has 

= x»D^ + ^(D^)4, (206) 

In order to obtain a covariant result, we have added in the last line two terms with 
Christoffel symbols whose sum cancels, as follows from 

= V + r F Y 

D$ = d^+'-u^^-h:^^ (207) 

We can, of course, add further covariant terms. A particular combination which 
will play a role is 

SsymW*) = <WX M ) +$eL(D[ m Xn]) 

D [m Xn] = \[{D mXn ) - (D nXm )\ (208) 

For the spinor field ip we then obtain 

1 ~ 1 \(D m y n ) - (D n y m )} 
S sym (X") ^ = X»D^ + -(D^W + ~ A - 1 2 1 n i m lni> (209) 
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To illustrate where this particular combination of Einstein and Lorentz transfor- 
mations comes from, we evaluate 5 sym on the vielbein. First we compute 5 cov on the 
vielbein 



&cov\X ^^m^ X ^u^m^ (fluXj^Tn X ^um 

= -dmX 11 + f(-r,p"e/) ^-D^ (210) 

where we defined d m = e m u d u and used the vielbein postulate. The transformation 
law defined in ( |208| ) then yields 

= -~(D mX » + D^xm) (211) 

Hence, if the vielbein would have been symmetric in its two indices m and fi to begin 
with, then a transformation with 5 sym keeps that symmetry. For this reason one 
might call the transformation 5 sym a "symmetric Einstein transformation", namely 
an Einstein transformation with parameter x^ which, as far as it acts on the vielbein, 
has been made symmetric by suitable local Lorentz transformations whose parameter 
also depends on x^- F° r spin-1/2 fields, both 5 cov in fl2U6| ) and 8 sym in fl209p are 



manifestly covariant, and at this point it is not clear which one to choose for the 
computation of gravitational anomalies. As we shall see, 5aw is actually equivalent 
to both. For spin-3/2, only one of them is covariant, as we now discuss. 

To understand the spin-3/2 transformation rule in (EU3), we first compute 5 cov 
on a spin-3/2 field ip m . We find 

= x^D^m + ^D^x^m (212) 

where again in the last line the two Christoffel symbols cancel. Hence, the proposed 
spin-3/2 law can be rewritten as 



^AW^m= &cov{X^) + 28iL [ D [mXn]) V>m (213) 



where the superscript (1) indicates that this Lorentz transformation only acts on the 
spin-1 (vector) index of if) m , but not on its spinor index. Compare this now with the 
symmetric spin-1/2 rule 

<W = [<Wx M ) + StL {1/2) (D [mX n])} $ (214) 
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where the superscript (|) indicates that this Lorentz transformation acts on the spinor 
indices. These rules for the spin-1/2 and spin-3/2 fields do not seem to agree at 
all. However, they are nevertheless equivalent, due to a surprising identity found by 
Fujikawa, Toniya, Yasuda and Endo [|17], [41| . Namely, if one uses pp as regulator 
both for the spin-1/2 case and for the spin-3/2 case, the anomalies due to |<Wj(x m ) 
are equal to minus those coming from Sil^ 2 ' {D[ m Xri\) '■ 

l^Uxn ~ -5eL {1/2 \D [mXn] ) (215) 
Using this equivalence, the proposed transformation rules for the spin-1/2 and the 



spin-3/2 field in ( [204]) turn out to be the same combination of Einstein and local 
Lorentz transformations, after all 

spin 1/2 : 5 AW = 6 cov = 25 cov + 25 eL 1/2 (Dx) = 25 sym 

spin 3/2 :<W = 5 cov + 25 lL {1 \D X ) 

= 25 cov + 26eL 1/2 (D X ) + 26 eL ^(Dx) 

= 25 sym (216) 

Note the relative factor 2 between 5 cov and 6 sym JKJ. 

It remains to prove that the anomalies coming from ^5 cov (x) are the same as 
those from —5il^(Dx)- We begin with the spin-1/2 fields. We take ip and ip m to 
be left-handed, whereas the fields i/j and tji m are independent right-handed fields in 
Euclidean space. (So ip is not ij}* in Euclidean space. In Minkowski space ip = ^'7°). 
We expand the spin-1/2 fields as follows: 

i> = °™ ^n,L + C a Xa,L ! $ = Y h n + Y d pCp,R j (217) 

n a n ft 

where ip n = g 1 ^ 4 ^ n and ip n are an orthonormal set of solutions of the Dirac equation 
with non- vanishing eigenvalue 

$J{x)$ n (x)dx = 5 mn , ipi> n = X n i> n , p = g 1/4 pg~ 1/4 (218) 

Of course, the chiral anomaly is equal to the number of minus the number of £'s, 
but we will rewrite this expression in a way where the difference between ip, x and ( 
disappears. 

For every eigenfunction ip n of the Dirac equation with nonvanishing eigenvalue 
X n , there is another one, 75^, with eigenvalue — A n . Hence, the "massive modes" 
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come in pairs, which can be combined on a chiral basis into ip n ,L and ip n ,Ri and 

i>mL^nRdx = , / IpmL^nLdx = / 1p m R^nRdx = 5 mn (219) 



where we used the orthogonality of i[) m and 7sV>n- (They are orthogonal because they 
are eigenfunction with different eigenvalues). In addition there are some zero modes, 
solutions of pip = 0- These we can always take to be left- and/or right-handed. They 
also appear in the expansion of ip and ?p, respectively, but they need not come in 
pairs. We have denoted them by Xa,L and C/3,_r in ( P17| )- 

The Jacobian J for an infinitesimal transformation 5ip and Sip follows from 



5a n = J dxipn^Sip , 5b n = J dx5ipip n , R (220) 
and similarly for 5c a and 8dp. For 5 aw we obtain then the following Jacobian 



Jaw 



dx E</W X*Dp+ -(D^ 



R 



XaL 
*Pn,R 



,3 



(221) 



where ip D^= d^ift — |'0^< < V m 7m7n,- The minus sign always appears in Jacobians 
for fermions. We now rewrite the massive modes in terms of ip n and 75^, but the 
zero modes we keep as they appear, except that we add a factor 75 or —75 to the 
left-handed or right-handed zero modes, respectively. (These factors 75 or —75 equal 
unity, of course). This yields 



Jaw~1= jdx^^n [x^+B^ + p^ 



4>n 



+ E ^75 



X »D,- D, X^ 



A. 



Jdx\y,x aL ] [x^D. + ^D^ 



IbXaL 



EW D,x" + ^D, X n 



(222) 



Partially integrating, all massive modes without factor 75 cancel, while the rest yields 



Jaw — 1 



/ dx 



+ Ea x« L h 5 [x"D^ + \ {D^)\ XaL + Ep <W7s [x^ + \{PpX?)\ (PR 
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(223) 



We would like to recognize this as a sum over a complete set of states. At first sight 
there is a mismatch in the first term which seems a factor 2 too large with respect to 
the other terms. However, the complete set of massive solutions contains not only ip n 
but also 75^n> an d rewriting 2ip n 'j 5 ip n for the massive modes as Vvys^n + i'-nlbi'-n 
with ip_ n = 75?/> n by definition, we arrive at 



Jaw - 1 = - J dx^, <Pn^5 



N 



2 



Vn (224) 



where (p^ is the complete set of eigenfunctions, so that N runs over all eigenf unctions: 
n > 0, n < 0, and all zero modes. 

To regulate the infinite sum over N, we add to each term a factor exp(— \ 2 N /M 2 ) 
and will let M tend to infinity at the end. Using 

e-^ViV = e ^' M2 VN (225) 



we obtain 



J AW - 1 = - f dxYs&N^n X^D^ + \{D^ 

N 1 



M 2 ~ 



(p N (226) 



where we recall the definition p = g x ^Jf)g 1 / 4 , with Jj) without Y^ y v term. All terms 
with r ^ cancel in ( |226|) , so from now on no explicit V ^ are present. 

Since the set ip^ is a complete set J2n l fN(x)<fiN^(y) = S(x — y), we have 



Jaw -I = -Tr ^-(x^ + D^e 



D„ = d^ + -^™ n lmln (227) 

where the trace Tr indicates integrating over space and summing over spinor indices. 

This trace has been evaluated using plane waves. Due to the g ±x ^ in Ip , the 
Einstein anomaly indeed did cancel [0, see below (|103| ). However, by using the 
cyclicity of the trace we can rewrite this as 

Jaw~1 = -Tr(i 75 (x^- 1/4 ^ 1/4 + g- 1/4 D„g^)e^ M2 ) (228) 
In the corresponding non-linear sigma model this becomes 

J AW - 1 = -L-Tr^xV* + 7T lX l ) exp(-(3H/h)) (229) 
Zirn 

with H given in ( |113| ) and 7Tj = — Y UJ iab' t P a " l P b ■ Note that the Jacobian is Weyl- 
ordered. One can from here compute the gravitational anomalies. 
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For local Lorentz transformations 

i 



1 



(230) 



one may proceed in a similar manner. One finds then for chiral spinors ip the following 
Lorentz anomaly 



JtL-l 



-Tr 7s 



4 A ' 



7m7rJ 



(231) 



Using cyclicity of the trace, all factors of g 1//4 and g 1//4 now cancel, and one obtains 

,1 



JtL-l = -Tr(- l5 \ mn lm7n exp(-(3H/h)) 



(232) 



Now that we have found expressions for the Einstein and Lorentz anomalies, we 
can prove the relation 5 cov ~ — \&tL for the spin- 1/2 fields. We use the identity 

/ Tr ( 75 (> + Pi)e^ /M2 ) = (233) 

which follows from cyclicity of the trace and ^75 = —j^p. Since the difference 
between pp as regulator and pp involves terms which are proportional to M~ 2 
in the path integral, and only the M-independent terms yield the anomaly while 
there are no singular terms in M for chiral anomalies, we may replace pp by pp 
in Q233|) . Then we can indeed use the cyclicity of the trace. (Chiral anomalies are 
rather insensitive to the details of the regulator. For trace anomalies, such details do 
matter, but for trace anomalies fortunately we do not need fl233|) ). Hence 



Tr 



Tr 
-2A 



dx 



dx 



dx 



(234) 



7 5 + D,xn + 757 V ] (x*Du - 

75 (X^D, + D^) + 757 ^(DpXv) 

mXn] ) 

This concludes the proof of the identity in ( |215| ) for spin- 1/2. 

For spin-3/2 we use the same identity in ([233), and all steps in (|234| ) are the 
same. Note that the local Lorentz transformations with D\ m Xn] only act on the 
spin-1/2 indices of the spin-3/2 field. Thus ( |215| ) is also proven for spin-3/2. The 
regulator pp now acts in a combined spinor-vector space, but it is diagonal in the 
vector indices, i.e., it acts on ip m the same for all m. 

To see whether this regulator is obtained from the gauge-fixed spin-3/2 action, 
we note that the spin-3/2 action in d dimension reads 



(235) 
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After adding a gauge fixing term [43 



£&< = d —^^^YYD^ p (236) 



and choosing a new basis for the spin-3/2 fields 



\fj, ~ ^7^7 ■ V> (237) 



the action becomes a sum of Dirac actions |42| 



£ + £ax = XnPx" = XmPx m (238) 

The field Xm transforms of course in the same way as ip m under space-time transfor- 
mations and the regulator for the spin-3/2 field \ m is thus pp, for the same reasons 
as for the spin-1/2 field. 



A. 4 Chiral anomalies in gravitational couplings from Feyn- 
man diagrams in quantum mechanics 

In this appendix we will derive the chiral anomaly in gravitational couplings using 
Feynman diagrams and the Feynman rules developed in the text. Starting point will 
be equations ( |126| ) and ( |129| ) derived in section 3.1. As explained there, we only have 
to consider one-loop diagrams. An arbitrary connected one-loop diagram consists of 
k vertices from (|129|) and k propagators. There are in principal many different ways 
to contract the vertices with each other, but using partial integration (in this case 
one can verify it is allowed) and the anti-symmetry of Rij a b in the indices i,j, one 
can always move the time derivative in (|129|) from qi to q % . This implies that all 
different contractions of the vertices yield the same contribution. The total number 
of contractions is 2 k ~ 1 (k — 1)!, which combines with the factor l/k\ from expanding 
( |129|) in the total symmetry factor 2 k ~ 1 /k for each diagram. If we always contract a 
q from one vertex with a q from another vertex, we find the following expression for 
5i ops, the set of connected one- loop diagrams 



1 oo nk—1 / ff\ k 

UrJ (239) 

where R k is R^Rll ■ ■ ■ R^ k with = g lk (xo)Rkj a b{u{x ))il)£ g ij)l g , and denotes the 
integral 



d n ... dr k (r 2 + 0(n - r 2 )) ■ • ■ (r k + 6^ - r fe ))(r a + 6(r k - n )) (240) 
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This expression is cyclically symmetric. To proceed, we compute YlkLi \ a k- This 
can be rewritten as 



y 



+ 



1 da 
o a 



a 



da 



^TL dn ---L dTk 



x(ar 2 + 6(n - t 2 )) ■ ■ ■ {ar k + 9{r k ^ - T k )){an + 6{r k - n))) . (241) 



Here, we introduced a parameter a in the expression for a k , which will turn out to 
make things much simpler. Notice that if we put this parameter equal to zero, we are 
left with an integral over a product of theta functions, which vanishes, except when 
k — 1, and then it equals 1/2 in view of 0(0) = 1/2, and this explains the first term 
y/2 in ( 241|) . The derivative d/da in ([241] ) can act on any of the k factors in the 
integral, but by cyclic symmetry these all give the same contribution and we choose 
it to hit only the last factor and put an extra factor of k in front. This yields 

v r 1 da ( 00 , r° r° \ 

f + 1 ^ ( E V J_ x dT l ■ ■ ■ j_ x ^(anXora + 0(n - r 2 )) • ■ • (ar k + 9{n-i - r k ))J . 

(242) 

Note that we now have broken the cyclic symmetry of ( |240| ) . If we expand the brackets 
inside the integral, we get a large sum of terms, each of which factors into a product 
of integrals of the form 



dr„ 



dr m+ i(aT m )9(T m - r m+i ) ■ • • #(r m+ ;_i - r m+i ). 



(243) 



For I = 0, the integrand is just ar m . This integral depends only on /, enabling us to 
write Will) as 




00 h f° f° 

Y.V dn... rfr fe (ari)0(r 1 - r 2 ) ■ • ■ 0(r fc _i - r k ) 
k=i J ~ l - 7 - 1 



(244) 



The integral in ( |243| ) equals —a/ {I + 2)!, and we find 



k=i K 




-a 



—a 



e y - 1 



l + s( e »_i_ y ) 



l0g (sin%/2)) 



This shows that 



^ 'S'loops 



- log 

2 6 



' TiR/A \ 
smh{TiR/A) ) 



(245) 
(246) 
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Thus, after a rescaling of the fermions, we find that the anomaly is given by 



An= Y[dx % y/g(x ) / n^bg ex P 



0=1 



1 / -iR/8n 

2 ° g Uinh(-iR/87r) 



(247) 



Clearly, there is only a gravitational 75 anomaly in n = 4k dimensions. In principle, 
this expression is ambiguous, since it is not clear whether R 4 means tr(i? 2 )tr(i? 2 ) or 
tr(i? 4 ). However, the precise meaning follows from the derivation above. First one 
has to write down a Taylor series for the logarithm, then one has to replace R m by 
tr(R m ) everywhere, and only then should one take the exponential. The advantage 
of writing the anomaly in the form ( 247|) is that it is given directly in terms of R, 
rather than in terms of its eigenvalues @ . 
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